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Abstract.  New  theoretical  methods  are  reported  for  ab  initio  calculations  of  the  adiabatic 
(Born-Oppenheimer)  electronic  wave  functions  and  potential  energy  surfaces  of  molecules 
and  other  atomic  aggregates.  An  outer  product  of  complete  sets  of  atomic  eigenstates 
familiar  from  perturbation-theoretical  treatments  of  long-range  interactions  is  employed 
as  a  representational  basis  without  prior  enforcement  of  aggregate  wave  function  anti¬ 
symmetry.  The  nature  and  attributes  of  this  atomic  spectral-product  basis  are  indicated, 
completeness  proofs  for  representation  of  antisymmetric  states  provided,  convergence  of 
Schrodinger  eigenstates  in  the  basis  established,  and  strategies  for  computational  imple- 
mention  of  the  theory  described.  A  diabatic-like  Hamiltonian  matrix  representative  is 
obtained  which  is  additive  in  atomic-energy  and  pairwise-atomic  interaction-energy  matri¬ 
ces,  providing  a  basis  for  molecular  calculations  in  terms  of  the  (Coulombic)  interactions 
of  the  atomic  constituents.  The  spectral-product  basis  is  shown  to  contain  the  totally 
antisymmetric  irreducible  representation  of  the  symmetric  group  of  aggregate  electron  co¬ 
ordinate  permutations  once  and  only  once,  but  to  also  span  other  (non-Pauli)  symmetric 
group  representations  known  to  contain  unphysical  discrete  states  and  associated  continua 
in  which  the  physically  significant  Schrodinger  eigenstates  are  generally  embedded.  These 
unphysical  representations  are  avoided  by  isolating  the  physical  block  of  the  Hamiltonian 
matrix  with  a  unitary  transformation  obtained  from  the  metric  matrix  of  the  explicitly  an¬ 
tisymmetrized  spectral-product  basis.  A  formal  proof  of  convergence  is  given  in  the  limit 
of  spectral  closure  to  wave  functions  and  energy  surfaces  obtained  employing  conventional 
prior  antisymmetrization,  but  determined  without  repeated  calculations  of  Hamiltonian 
matrix  elements  as  integrals  over  explicitly  antisymmetric  aggregate  basis  states.  Com¬ 
putational  implementations  of  the  theory  employ  efficient  recursive  methods  which  avoid 
explicit  construction  the  metric  matrix  and  do  not  require  storage  of  the  full  Hamilto¬ 
nian  matrix  to  isolate  the  antisymmetric  subspace  of  the  spectral-product  representation. 
Calculations  of  the  lowest-lying  singlet  and  triplet  electronic  states  of  the  covalent  elec¬ 
tron  pair  bond  (H2)  illustrate  the  various  theorems  devised  and  demonstrate  the  degree  of 
convergence  achieved  to  values  obtained  employing  conventional  prior  antisymmetrization. 
Concluding  remarks  place  the  atomic  spectral-product  development  in  the  context  of  cur¬ 
rently  employed  approaches  for  ab  initio  construction  of  adiabatic  electronic  eigenfunctions 
and  potential  energy  surfaces,  provide  comparisons  with  earlier  related  approaches,  and 
indicate  prospects  for  more  general  applications  of  the  method. 
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1.0  Introduction 


Adiabatic  electronic  wave  functions  provide  well-known  first  approximations  to  the 
physical  attributes  of  molecules  and  other  forms  of  matter,1-3  and  are  prerequisites  for 
studies  of  the  non-adiabatic  dynamical  consequences  of  avoided  crossings  and  conical  inter¬ 
sections  in  the  potential  energy  surfaces  that  guide  the  pathways  of  chemical  reactions.4-6 
Methods  currently  employed  for  constructing  such  wave  functions  and  energy  surfaces 
include  ab  initio  quantum  chemistry,7-10  quantum  Monte-Carlo  simulations,11’12  density- 
functional  approaches,13-18  semi-empirical  techniques,19-23  hybrid  quantum  mechanical/ 
molecular  mechanical  (QM/MM)  combinations,24  and  various  forms  of  long-range  pertur¬ 
bation  theory.25-38  Although  these  methods  have  provided  a  significant  amount  of  useful 
information  over  the  past  few  decades,  they  arguably  do  not  constitute  a  quantum  theory 
of  chemical  bonding,  which  would  preferably  express  the  chemical  attributes  of  compounds 
in  terms  of  intrinsic  properties  of  their  atomic  constituents.  New  approaches  in  this  spirit 
to  construction  of  adiabatic  electronic  wave  functions  are  clearly  welcome,  particularly  if 
they  can  retain  the  desirable  features  of  the  existing  methods  while  ameliorating  some  of 
their  shortcomings. 

The  ab  initio  approaches  employing  multi-configurational  Hartree-Fock,  Mpller-Plessett 
perturbation  theory,  and  multi-reference  configuration  mixing  methods,7,8  and  the  coupled- 
cluster,9  valence-bond,10  and  quantum  Monte  Carlo11’12  approaches,  can  provide  highly 
accurate  energies  and  other  properties  for  the  ground  and  excited  states  of  small  molecules, 
but  they  generally  require  allocations  of  significant  computational  resources  for  this  pur¬ 
pose,  and  they  are  apparently  not  yet  applicable  to  large  atomic  aggregates.  By  con¬ 
trast,  the  widely  employed  density-functional,13-18  semi-empirical,19-23  and  QM/MM24 
approaches  can  be  applied  to  larger  aggregates,  including  atomic  clusters,22  condensed 
matter  systems,23  and  biological  macromolecules,24  and  in  certain  cases  provide  electron¬ 
ically  excited  states,  although  their  a  priori  accuracy  is  frequently  unknown.  Long-range 
Rayleigh- Schrodinger  and  symmetry- adapted  perturbation  theories  have  also  proved  useful 
in  appropriate  limits,25-38  although  widespread  and  systematic  applications  of  perturba¬ 
tion  theory  to  chemical  interactions  are  apparently  not  yet  in  evidence.  These  methods 
have  all  contributed  to  the  collective  development  of  a  significant  methological  and  compu¬ 
tational  base  from  which  to  draw  upon  in  devising  new  ab  initio  approaches  to  electronic- 
structure  and  spectral  calculations.39-44 

In  the  present  manuscript,  new  theoretical  methods  are  described  for  ab  initio  calcula¬ 
tions  of  the  adiabatic  electronic  wave  functions  and  potential  energy  surfaces  of  molecules 
and  other  atomic  aggregates.45  Certain  of  the  desirable  attributes  of  the  computational 
methods  currently  in  play  for  these  purposes  are  incorporated  in  the  development,  which 
employs  an  outer  product  of  complete  sets  of  atomic  spectral  eigenfunctions  familiar  from 
the  aforementioned  perturbation-theoretical  treatments  of  long-range  atomic  interactions.25 
In  contrast  to  related  earlier  developments  which  employ  prior  antisymmetrization  of  prod¬ 
ucts  of  atomic  wave  functions  in  atoms-in-molecules  approaches,46-58  antisymmetry  is 
enforced  in  the  present  development  subsequent  to  construction  of  the  Hamiltonian  ma¬ 
trix  employing  a  unitary  transformation  to  isolate  its  antisymmetric  subspace.  The  to¬ 
tal  electronic  Hamiltonian  operator  in  the  atomic-product  representation  takes  a  matrix 
form  in  which  the  individual  atomic  energies  appear  explicitly  as  constants  independent 
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of  atomic  positions,  with  only  a  sum  of  pairwise- atomic  interaction-energy  matrices  re¬ 
quiring  repeated  evaluation  as  a  function  of  aggregate  atomic  configuration  in  potential- 
energy- surf  ace  calculations.  The  pairwise  interaction-energy  matrices  depend  only  upon 
the  vector  separations  of  the  two  interacting  atoms,  are  independent  of  the  positions  of  all 
other  atoms  in  the  aggregate,  and  can  be  expressed  in  terms  of  one-electron  atomic  charge 
and  transition  densities,  which  provide  the  computational  invariants  of  the  development. 
Accordingly,  the  present  approach  does  not  require  repeated  calculations  of  Hamiltonian 
matrix  elements  as  integrals  over  antisymmetric  many-electron  aggregate  basis  functions; 
rather,  one-electron  atomic  charge  and  transition  densities  are  calculated  once  and  for 
all  for  a  spectrum  of  atomic  states  and  retained  for  repeated  use  in  construction  of  the 
Hamiltonian  matrix  for  arbitrary  spatial  arrangements  of  interacting  atoms. 

Notational  conventions  are  established  and  the  (Coulombic)  form  of  the  many-electron 
Hamiltonian  operator  employed  in  the  development  given  in  Section  2.  The  atomic  spectral- 
product  basis  and  certain  of  its  attributes  are  described  in  Section  3  employing  complete 
sets  of  spin  orbitals  in  representations  of  the  many-electron  atomic  spectral  eigenstates.  It 
is  shown  that  the  totally  antisymmetric  irreducible  representation  of  the  aggregate  symmet¬ 
ric  group  is  spanned  once  and  only  once  by  the  atomic-product  basis,  although  other  (non- 
Pauli)  symmetric-group  representations  are  also  contained  in  the  reducible  product.  The 
form  of  the  Hamiltonian  matrix  in  the  spectral-product  representation  is  given  in  Section  4, 
methods  are  described  for  its  evaluation  in  terms  of  one-electron  atomic  transition-density 
matrices,  and  certain  attributes  of  the  spectrum  of  the  Hamiltonian  matrix  are  indicated. 
In  Section  5,  methods  are  reported  for  formal  isolation  of  the  totally  antisymmetric  block  of 
the  spectral-product  Hamiltonian  matrix  employing  a  unitary  transformation  derived  from 
the  metric  matrix  of  the  antisymmetrized  spectral-product  basis.  Convergence  is  demon¬ 
strated  in  this  way  to  eigenstates  and  energies  identical  to  those  obtained  in  the  linearly 
independent  subspace  of  this  prior  antisymmetrized  representation.  Aspects  of  the  theory 
are  illustrated  with  calculations  reported  in  Section  6,  which  employ  efficient  recursive 
methods  for  isolating  the  totally  antisymmetric  block  of  the  Hamiltonian  matrix,  avoiding 
explicit  construction  of  the  metric  matrix  and  of  the  entire  Hamiltonian  matrix  at  any  one 
time.  Convergence  to  the  lowest-lying  singlet  and  triplet  potential  energy  curves  of  the 
electron  pair  bond  (H2)  is  achieved  employing  even-tempered  Gaussian  orbital  representa¬ 
tions  and  related  computational  methodology.  Finally,  the  new  development  is  discussed  in 
the  context  of  currently  employed  electronic-structure  approaches,  comparisons  are  made 
with  related  earlier  theoretical  methods,  and  prospects  for  more  general  applications  of 
the  method  are  indicated  in  concluding  remarks  provided  in  Section  7. 

2.0  The  Aggregate  Schrodinger  Equation 

The  adiabatic  electronic  eigenstates  and  energies  of  an  iV-atom  molecule  or  other  atomic 
aggregate  are  obtained  from  the  Schrodinger  equation 


F7(i,  2, . . . ,  n  :  R)^( l,  a, . . . ,  n  :  R)  —  ^(l,  2, . . . ,  n  :  R )  •  E (R), 


(1) 


where  \l/(i,  2 , ,n  :  R)  is  a  row  vector  of  orthonormal  eigenstates  \lrr(i,2,...,n  :  R) 
having  quantum  labels  T  and  E (R)  is  the  diagonal  matrix  of  corresponding  eigenvalues 
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E-p(R)-  The  Hamiltonian  operator  appearing  in  Eq.  (1)  is  written  in  the  form 


N  N 

£(i,  a,..,,  n:K)  =  £{£<“>(<  :fl„)+  £  j  :  R^)}, 


(2) 


a=l 


/3=a  +  l 


where 
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H^(i:Ra)  =  J2{ 
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na 

+  £ 

*'=*+! 


(3) 


is  the  Hamiltonian  operator  in  the  familar  Coulomb  approximation  for  a  set  of  na  electrons 
which  are  arbitrarily  associated  with  the  atomic  nucleus  at  position  Rai  and  whose  3na 
spatial  coordinates  are  represented  by  the  single  vector  label  i,  and 
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is  the  Coulombic  interaction  potential  between  atoms  a  and  ft.  where  the  coordinate  label  j 
refers  to  the  3 np  spatial  electron  coordinates  arbitrarily  associated  with  the  atomic  nucleus 
at  the  position  Rp.59  The  indecies  l,  2, . . . ,  n  in  Eqs.  (1)  and  (2)  represent  the  collected 
spatial  coordinates  of  disjoint  sets  of  electrons,  each  arbitrarily  associated  with  one  of  the 
atomic  nuclei  in  the  aggregate. 

In  the  foregoing  equations  and  throughout  the  development,  all  electron  (rft  and 
atomic-position  (Ra)  coordinates  are  defined  in  the  laboratory  frame,  rja  =  r*  —  Ra 
and  r:lfj  =  r:,  —  Rp  are  electron  coordinate  vectors  relative  to  the  indicated  atomic  po¬ 
sitions,  rij  —  ri  —  rj  is  the  (i.  j)  electron  vector  separation,  Rap  =  Ra  —  Rp  is  the 
(a,  ft)  atomic  position  vector  separation,  R  =  (Ri,  R2, . . . ,  Rn)  specifies  the  positions 
of  all  the  atoms  in  the  entire  aggregate,  e  is  the  magnitude  of  the  electronic  charge,  m 
is  the  electron  mass,  Z(x  and  Zp  are  atomic  numbers,  and  the  other  symbols  have  their 
usual  meanings.59  A  colon  is  employed  to  separate  electron  coordinates  from  atomic  po¬ 
sitions  (i  :  Ra ),  and  a  semi-colon  is  employed  in  Eq.  (4)  to  distinguish  sets  of  electrons 
(i:  j)  between  which  the  interaction  potential  :  Rap)  is  not  totally  symmetric 

under  coordinate  permutations.  Of  course,  the  n^-electron  (nt  =  n\  +  n2  +  ■  ■  ■  tin)  aggre¬ 
gate  Hamiltonian  H(i,  2, . . .,  n  :  R)  and  the  individual  na-electron  atomic  Hamiltonians 
H^{i  :  Ra)  are  symmetric  sums  of  one-  and  two-electron  operators  in  their  respective 
electron  labels.  Accordingly,  the  notation  employed  in  Eqs.  (1)  to  (4)  does  not  by  itself 
constitute  a  meaningful  assignment  of  particular  electrons  to  particular  atoms,  although 
it  does  suggest  this  possibility  in  conjunction  with  the  choice  of  an  appropriate  electronic 
representational  basis. 

The  physical  eigenstates  \lAp)(i,  2, . . . ,  n  :  R)  and  eigenvalues  E ^(R)  included  in  the 
solutions  of  Eq.  (1)  for  the  Hamiltonian  of  Eqs.  (2)  to  (4)  have  good  total  spin  and  spin- 
projection  quantum  numbers  (S,  Ms)  and  transform  as  the  totally  antisymmetric  repre¬ 
sentation  (lnt)  of  the  symmetric  group  Snt  of  the  aggregate  electron  spin  and  spatial  coor¬ 
dinate  permutations,60  whereas  additional  (non-Pauli)  eigensolutions  2, . . . ,  n  :  R) 
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of  Eq.  (1)  which  transform  as  other  than  the  totally  antisymmetric  irreducible  representa¬ 
tion  under  electron  coordinate  permutations  are  generally  regarded  as  unphysical,  as  are 
the  corresponding  energies  E^)(i?).61 

3.0  Definition  and  Attributes  of  the  Spectral-Product  Basis 

The  physical  and  unphysical  solutions  of  Eq.  (1)  which  together  comprise  the  full 
spectrum  of  the  Hamiltonian  operator  of  Eqs.  (2)  to  (4)  can  be  constructed  in  a  many- 
electron  spectral-product  basis  the  nature  and  attributes  of  which  are  described  in  this 
section.  The  presence  of  unphysical  states  and  the  importance  of  devising  methods  for 
their  elimination  is  emphasized  throughout  the  development. 

3.1  Spin-Orbital  Basis  Sets 

Conventional  variational  approximations  to  the  solutions  of  Eq.  (1)  in  their  most 
common  form  employ  one-electron  basis  sets  of  atomic  spin  orbitals  in  Slater,  Gaussian, 
or  other  forms  centered  at  each  of  the  aggregate  atoms  (a  —  1,2,...,  N).  Such  basis  sets, 
employed  here  largely  for  purposes  of  analysis,62  are  conveniently  written  as  row  vectors 
(f>(a)(i  :  Ra  )  of  orthonormal  one-electron  functions  (i  :  Ra)  of  the  individual  electronic 
spin  and  of  spatial  coordinates  measured  relative  to  the  atomic  positions  (i  =  r?;  —  Ra). 
These  functions  are  generally  enumerated  employing  an  appropriate  set  of  quantum  labels 
7a,  are  assumed  complete  in  a  finite  three-dimensional  spatial  domain,  and,  accordingly, 
satisfy  the  closure  relation 


OO 

0(o)(* :  Ra)  '  0<a)(»'  :  flat  =  £  4>^(i  :  <* (a)(*  -  *')  (5) 

7«  =  1 

in  this  domain,  where  :  Ra )t  is  a  column  vector  comprised  of  the  complex  conjugate 

spin  orbitals  (f>^J (i'  :  Ra)* .  The  choice  of  particular  spin  orbitals  employed  in  Eq.  (5) 
is  arbitrary  in  the  limit  of  a  complete  set,  but  the  rate  of  convergence  of  computational 
developments  built  upon  them  will  generally  depend  upon  the  specific  basis  employed.62 

Traditional  uses  of  these  one-electron  basis  sets  in  connection  with  solution  of  Eq. 
(1)  follow  well-known  molecular-orbital  or  valence-bond  methods  in  constructing  many- 
electron  representational  basis  sets.63  These  conventional  developments  have  in  common 
the  representation  of  each  electron  in  a  “symmetrical”  or  “democratic”  manner  and  the  use 
of  prior  term-by-term  electron  coordinate  antisymmetrization  of  the  many-electron  basis 
employed.61  In  the  molecular- orbital  approaches,  the  common  one-electron  basis  set  that 
extends  over  the  3rq-dimensional  spatial  domain  of  the  aggregate  is  employed  to  describe 
the  positions  of  any  of  the  electrons  in  antisymmetric  Slater  determinants,7’8  whereas  in  the 
valence-bond  methods  over-all  term-by-term  prior  antisymmetry  insures  the  canonical  spin- 
paired  product  states  employed  treat  the  electrons  in  a  democratic  fashion.10  Of  course, 
use  of  prior  antisymmetry,  and  of  symmetry  adaptation  more  generally,  guarantees  that 
the  matrix  representative  of  the  Hamiltonian  operator  obtained  will  be  block  diagonal  and 
of  minimal  possible  dimension  for  the  states  of  interest  in  the  basis  employed,  in  accordance 
with  Schur’s  lemmas.60 

3.2  Many-Electron  Outer-Product  Basis  Set 
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In  contrast  to  the  aforementioned  conventional  developments,  it  is  possible  to  forego 
prior  symmetry  adaptation  of  the  individual  terms  in  the  many-electron  basis  set  in  favor 
of  other  potential  advantages  and  simplifications  by  employing  ^-electron  states  which 
are  not  individually  antisymmetric  at  the  outset.  Such  states  can  be  constructed  from 
products  of  na-electron  atomic  basis  states  which  themselves  are  written  in  the  form  of 
simple  products  of  spin  orbitals, 


^  a 

7« 

where  the  label  Ta  specifies  a  particular  choice  of  spin  orbitals  on  the  atom  a  in  the  indi¬ 
cated  na-term  product.  The  row  vector  Il(a)(i  :  Ra)  obtained  from  all  Ta  values  includes 
all  possible  unrestricted  spin-orbital  products  in  Eq.  (6),  thereby  providing  a  complete 
representational  basis  for  na-electron  eigenstates  consequent  of  Eq.  (5).  Accordingly,  the 
complete  outer  product  of  such  states  for  all  atoms  in  the  aggregate,  written  in  the  form 
of  the  complete  separable  Hilbert  space 

II(i;  2; . . . ;  n  :  R )  =  jll^i  :  Rt)  <g)Il(2)(2  :  R2)  ®  ■  •  •  n^(n  :  RN)^  (7) 

similarly  provides  a  complete  representational  basis  for  the  entire  nrelectron  aggregate. 
Here,  semi-colons  are  employed  to  separate  groups  of  electrons  (i;  a; . . . ;  n)  distinguished 
by  virtue  of  their  representations  in  the  different  sets  of  spin-orbital  product  states  em¬ 
ployed  for  each  atom,  whereas  the  outer- product  symbol  ”<g>”  is  used  to  indicate  that  all 
possible  products  of  the  atomic  basis  functions  are  included  in  the  vector  II(i;  a; . .  471  : 
R).  The  bracket  symbol  {•  •  -}Q  implies  choice  of  a  particular  rule  for  enumerating  and 
ordering  the  product  functions  in  the  row  vector  II (1;  a; . . . ;  n  :  R ),  which  rule,  although 
arbitrary,  must  be  adhered  to  in  matrices  constructed  in  the  outer-product  basis  and  in 
related  vector  and  matrix  multiplications. 

The  basis  of  Eq.  (7),  in  which  specific  sets  of  electrons  (1;  a; . . . ;  n)  are  arbitrarily  as¬ 
signed  to  specific  atoms  (1,2,...,  N )  and  are  described  employing  different  sets  of  orbital 
functions  [Eq.  (5)],  does  not  provide  individual  rq-electron  functions  which  transform  irre- 
ducibly  under  aggregate  electron  coordinate  permutations,  nor  are  these  functions  generally 
eigenstates  of  total  electron  spin  or  of  other  aggregate  commuting  operators.  Nevertheless, 
as  a  consequence  of  the  closure  relations  of  Eq.  (5)  for  the  different  sets  of  spin  orbitals, 
the  aggregate  outer-product  basis  satisfies  the  closure  relation 


N 


N 


n(i;  2; . . . ;  n  :  R)  U( 1';  2'; . . . ;  n'  :  R)^  =  J]  Il(a)(i  :  Ra)-U^\i'  :  R^  ]J  5<a) (*-*'), 


a= 1 


a= 1 


(8) 

where  the  product  of  delta  functions  on  the  right-hand  side  expresses  completeness  of  the 
basis  as  a  separable  Hilbert  space  in  the  spin  coordinates  and  over  a  3nt-dimensional  spa¬ 
tial  domain.  Accordingly,  convergence  at  least  in  the  norm  to  physical  and  unphysical 
eigenstates  of  irreducible  symmetry  can  be  achieved  in  the  basis  of  Eqs.  (5)  to  (8),  poten¬ 
tially  providing  approximations  to  any  of  the  aggregate  Schrodinger  eigensolutions  of  Eq. 
(!)• 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


3.3  Atomic  Spectral  Eigenstates 

It  is  convenient  to  perform  a  series  of  formal  unitary  transformations  of  the  complete 
product  basis  sets  of  Eq.  (6)  which  bring  the  matrix  representatives  of  the  individual 
atomic  Hamiltonian  operators  of  Eq.  (3)  into  diagonal  forms,  thereby  simplifying  the 
total  Hamiltonian  matrix  but  leaving  the  solutions  of  Eq.  (1)  in  the  basis  invariant. 
Specifically,  because  a  single  spin-orbital  basis  set  is  used  for  all  electrons  i  on  atom  a, 
and  all  spin-orbital  products  are  present  in  Eq.  (6),  all  permutations  among  electron 
coordinates  i  are  also  included  explicitly  for  each  atomic  site.  It  should  be  recalled  in  this 
connection  that  all  permutations  of  electron  coordinates  among  the  spin-orbital  products  of 
Eq.  (6)  are  included  explicitly  in  conventional  atomic  structure  theory,  so  that  the  orbital 
configurations  employed  to  represent  atomic  states  must  be  limited  to  sets  of  ordered  spin- 
orbital  products  in  Slater  determinants  to  insure  the  linear  independence  of  the  basis  states 
employed.62  The  present  development,  in  which  all  orderings  of  spin-orbital  products  for 
a  given  atom  are  included  but  only  a  single  assignment  of  electrons  to  spin  orbitals  is 
made,  is  entirely  equivalent  to  the  standard  method  because  a  common  spin-orbital  basis 
is  employed  in  a  symmetrical  representation  of  all  the  electrons  in  a  given  atom.61 

The  eigenstates  of  the  atomic  Hamiltonian  of  Eq.  (3), 

H^a\i  :  Ra)^a\i  :  Ra)  =  :  Ra)  •  E^,  (9) 

are  represented  in  the  form 

&a\i:Ra)  =  Il(a\i:Ra)-U(a\  (10) 

where  the  unitary  transformation  matrix  Tj(a)  is  independent  of  Ra,  and  E(a)  is  the  diag¬ 
onal  matrix  of  all  atomic  eigenvalues  of  the  atom  a  obtained  from  the  matrix  Schrodinger 
equation 

<n(a)(i  :  Ra)\H^a\i  :  Ra)\U^a\i  :  Ra))  •  U<a>  =  U<a>  •  eK  (11) 

In  Eq.  (11),  an  outer-vector-product  convention  is  employed  on  the  left-hand  side  to  pro¬ 
vide  matrices  having  conventional  row  and  column  labels  formed  from  the  indicated  bra 
(n<a)(i  :  Ra)  |  and  ket  |Il(a)(i  :  Ra))  vectors.  The  irreducible-symmetry  labels  or  good 
quantum  numbers  of  the  eigenstates  of  Eq.  (9)  can  also  be  accommodated  into  the  trans¬ 
formation  matrix  of  Eqs.  (10)  and  (11)  without  loss  of  generality,  and  so  <!>(“)  (i  :  Ra)  can 
be  taken  to  be  a  row  vector  of  functions  ^[^(i  :  Ra)  which  are  multiplet  eigenstates  of 
the  Hamiltonian  H^a\i  :  Ra )  having  specific  energy  Era,  orbital  and  spin  angular  mo¬ 
mentum,  and  parity  values;  i.e.,  Ta  =  (E,T,M£,S,,Ms,P)a.63  These  atomic  eigenfunctions 
transform  either  as  physical  (totally  antisymmetric)  or  unphysical  (non-Pauli)  irreducible 
representations  under  the  coordinate  permutation  group  SUa  for  electrons  in  the  set  i.60 

Combining  the  unitary  transformations  of  Eqs.  (10)  and  (11)  for  all  atoms  in  the 
aggregate,  a  complete  set  of  products  of  atomic  spectral  eigenstates  is  obtained  in  the 
form 

<h(i;  2; . . . ;  n  :  R)  =  II(i;  2; . . . ;  n  :  R)  ■  U,  (12) 
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where  U  =  {U^  (8)  U^2)  (8)  •  •  •  U^}0  is  the  appropriately  ordered  outer  matrix  product 
of  unitary  transformation  matrices  of  the  individual  spectral  basis  sets  of  Eqs.  (10)  and 
(11),  and 

S(i;a;...;n:  R)  =  j$(1)(i  :  Ri)  <8>  #(2)(2  :  R2)  <8>  ■  ■  •  &N\n  :  ^iv)}Q  (13) 

is  the  resulting  complete  set  of  ordered  products  of  atomic  spectral  eigenstates.  As  a 
consequence  of  Eq.  (8)  and  the  definitions  of  Eqs.  (9)  to  (12),  the  spectral-product  states 
of  Eq.  (13)  satisfy  the  closure  relation 


N  N 

S(i;a;...;n:  R)-$(i';  a'; . . , ;  n'  :  i?)t  =  JJ  $(a)(i  :  Ra)-&a)(i'  ■  RaV  ->•  ]J  5(a)(*-*') 

a=l  a= 1 

(14) 

for  all  irreducible  representations  of  S„t,  and  so  can  be  employed  in  descriptions  of  any 
of  the  eigensolutions  of  Eq.  (1).  Use  of  the  spectral-product  representation  of  Eq. 
(13)  in  place  of  the  outer-product  representation  of  Eq.  (7)  constitutes  a  partial  pre- 
diagonalization  of  the  aggregate  Hamiltonian  matrix  which  leaves  its  dimension  unchanged 
and  otherwise  provides  results  equivalent  to  those  obtained  from  the  outer-product  basis. 

3.4  Restriction  to  Antisymmetric  Atomic  States 

In  view  of  the  non-symmetrical  nature  of  the  spectral-product  basis  of  Eq.  (13),  and 
of  the  fact  that  the  outer  product  of  the  symmetric  groups  of  the  individual  atomic  con¬ 
stituents  forms  a  subgroup  of  the  full  electron-coordinate  permutation  group  ( <Si <8>  S2  <8>  •  •  • 
Sn  C  Sth ) ,  some  care  must  be  exercised  in  employing  standard  group-theoretical  proce¬ 
dures  in  identifying  and  isolating  the  desired  totally  antisymmetric  representation  (lnt)  of 
S„t  contained  in  the  basis.  The  appropriate  group-theoretical  product  in  this  case  is  the 
outer  product  of  subgroup  representations,  as  opposed  to  the  inner  Kronecker  group  prod¬ 
uct  which  is  commonly  reduced  employing  a  Clebsch-Gordan  decomposition.60  Moreover, 
in  the  representation  of  Eq.  (13),  although  electron-coordinate  permutations  within  each 
set  i  are  included  explicitly ,  those  permutations  involving  at  least  one  transposition  of  elec¬ 
tron  coordinates  from  different  sets  i  and  j  at  the  two  atomic  sites  a  and  fd  are  spanned 
implicitly.  Nevertheless,  in  the  limit  of  the  closure  associated  with  use  of  all  orderings 
of  spin-orbital  products  in  Eq.  (13),  the  latter  implicit  electron  coordinate  permutations 
are  all  effectively  present  in  Eq.  (13)  and  the  standard  outer- product  reduction  theory 
applies.60 

Although  there  is  apparently  no  convenient  way  in  which  to  use  group-theoretical  meth¬ 
ods  in  the  explicit  construction  of  a  particular  irreducible  subspace  in  the  ordered-product 
representation  of  Eq.  (13),  the  outer-product  reduction  procedure  can,  nevertheless,  be 
employed  as  an  analytical  device  in  this  connection.32  In  particular,  it  can  be  demon¬ 
strated  that  some  (but  not  all)  of  the  unwanted  unphysical  representations  of  S„t  spanned 
by  the  basis  of  Eq.  (13)  are  removed  by  the  simple  expedient  of  restricting  the  atomic 
eigenstates  in  the  vectors  <lAa)(i  :  Ra)  for  all  atoms  a  to  the  physical  subspace  “p”  of 
eigenstates  :  Ra)  which  are  antisymmetric  in  the  na  electron  coordinates  i.  The 

required  demonstration  follows  from  the  standard  outer-product  reduction  rules,60  which 
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indicate  that  the  presence  of  any  horizontal  (“symmetric”)  strip  in  either  of  two  Young 
patterns  associated  with  distinct  symmetric  groups  disallows  the  appearance  of  the  totally 
antisymmetric  representation  in  the  reduction  of  the  outer- product  of  the  two  patterns  in 
the  combined  permutation  group.  More  specifically,  the  reduction  of  the  outer  product  of 
any  two  irreducible  representations  (7n  &jm)  of  two  distinct  symmetric  groups  Sn .  Srn  into 
irreducible  representations  of  the  group  Sn+m  includes  the  totally  antisymmetric  represen¬ 
tation  i(n+m)  of  Sn+m  once  and  only  once  if  and  only  if  the  two  subgroup  representations 
are  the  totally  antisymmetric  representations  7n  =  ln  and  7 m  =  lm.  Letting  Sn  represent 
the  symmetric  group  for  the  first  a  —  1  atoms  in  the  aggregate  (n  =  ni+ri2  +  -  •  •  na_  1),  set¬ 
ting  Sm=Sna,  the  symmetric  group  of  the  atom  o;,  and  assuming  that  1”  is  contained  only 
once  in  Sni  it  is  seen  that  ln  ®  r)Ua  contains  l(n+n<»)  only  once  if  and  only  if  /jna  =  ln“. 
Correspondingly,  it  follows  by  induction  that  the  totally  antisymmetric  representation 
lnt  of  S„t  is  obtained  once  and  only  once  in  the  reduction  of  the  entire  outer  product 
(7m  ®  7n2  ®  ■  ■  ■  7 nN)  if  and  only  if  ^Ua  =  ln“  for  all  a. 

In  accordance  with  the  foregoing  observations,  a  complete  but  not  over-complete  linear ly- 
independent  orthonormal  subspace  of  the  spectral-product  basis  of  Eq.  (13)  suitable  for 
representations  of  totally  antisymmetric  aggregate  states  is  contained  in  the  products  of 
the  physical  (antisymmetric)  atomic  eigenstates  <hpa)(i  :  Ra), 

$(1;  2; . . . ;  N  :  R)  =  {^(i  :  Ri)  ®  $<2)(2  :  R2)  ®  •  •  •  <f>{pN\n  :  i?iV)}o  •  (15) 

The  spectral  closure  relation  of  Eq.  (14)  correspondingly  now  becomes  operative  only 
in  the  subspace  of  irreducible  representations  of  Snt  spanned  by  the  outer  product  of  all 
antisymmetric  atomic  eigenstates, 


N 


N 


S(i;  =»j  ■ .  • ;  n  :  3';  =  ff  :  :  Ra)'  -+  ff 


a= 1  a=l 

(16) 

where  Spa\i  —  i')  refers  to  the  Dirac  delta  function  in  the  totally  antisymmetric  physical 
subspace  (i  :  Ra)  of  the  atomic  states  [cf.,  Eq.  (14)]. 

Although  some  of  the  non-totally-antisymmetric  representations  of  Snt  present  in  the 
outer-product  basis  are  generally  absent  in  the  basis  of  Eq.  (15),  a  certain  number  of 
these  non-Pauli  representations  are  always  present.  The  number  and  type  of  unphysical 
irreducible  representations  of  Snt  spanned  by  the  reducible  spectral-product  representation 
of  Eq.  (15)  is  determined  by  the  total  number  of  electrons  nt  and  by  the  partitioning 
of  these  into  the  particular  atomic  values  n\,  ri2,  ■  ■  ■ ,  un  employed.  Although,  strictly 
speaking,  these  unphysical  irreducible  representations  are  not  explicitly  required  in  the 
present  development,  they  can  be  enumerated  in  specific  cases  employing  the  outer- product 
reduction  rules  in  the  usual  way.60  By  this  identification  it  is  made  clear  that  the  basis 
of  Eq.  (15)  can  span  unphysical  states  of  very  low  energy  arising  from  “over- filled”  inner 
shells  and  their  associated  continua,  in  which  the  physical  states  can  be  embedded.32-37 

The  basis  of  Eq.  (15)  has  been  employed  previously  in  a  number  of  related  connec¬ 
tions,  including  long-range  and  exchange  perturbation  theories25-38  and  so-called  atoms- 
in- molecules46- 50  and  related  diatomics-in-molecules51-58  methods.  In  contrast  to  the 
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present  development,  the  atoms-in-molecules  approaches  employ  explicit  prior  antisym- 
metrization  of  the  spectral-product  basis  of  Eq.  (15),  in  spite  of  the  consequent  linear 
dependence  in  the  limit  of  closure  this  approach  entails,  and  the  perturbation-theoretical 
methods  are  generally  applicable  only  in  cases  of  weak  long-range  atomic  interactions.  The 
representation  of  Eq.  (15)  would  appear  to  be  particularly  appropriate  for  molecules  and 
other  atomic  aggregates  even  in  the  absence  of  prior  aggregate  antisymmetry  in  that  it 
is  strictly  orthonormal  and  contains  the  antisymmetric  representation  only  once,  avoiding 
thereby  problems  of  linear  dependence.  It  furthermore  isolates  the  most  singular  intra- 
atomic  Coulombic  interactions  in  the  individual  atomic  wave  functions  employed  as  a  basis, 
gives  the  exact  aggregate  energies  in  complete  dissociation  limits,  and  results  in  a  Hamil¬ 
tonian  matrix  of  particularly  simple  form  (see  below).  Additionally,  the  basis  of  Eq.  (15) 
provides  the  possibility  of  developing  an  atomic-interaction-based  approach  for  evaluating 
an  aggregate  eigenspectrum  in  which  the  exclusion  principle  among  aggregate  electrons  is 
accommodated  in  accordance  with  the  degree  of  interaction  between  atoms,  rather  than  as 
a  global  constraint  applied  uniformly  and  uncritically  among  all  electrons  present,  which 
constraint  is  generally  not  required  in  the  case  of  highly  separated  or  otherwise  weakly- 
interacting  atomic  constitutents.64 

4.0  Spectral-Product  Representation  of  Schrodinger  Eigenstates 

The  spectral-product  basis  defined  in  the  preceding  Section  can  be  employed  to  con¬ 
struct  solutions  of  Eq.  (1)  following  standard  variational  procedures.  Specifically,  employ¬ 
ing  an  expansion  of  eigenfunctions  in  the  spectral-product  basis  of  Eq.  (15)  in  Eq.  (1)  and 
projecting  with  each  of  the  orthonormal  product  functions  in  the  basis  in  the  usual  way 
gives  the  familiar  matrix  Schrodinger  equation63 

H (R)  •  UH (R)  =  Uh(H)  •  E (R),  (17) 

where  E(i?)  is  the  desired  diagonal  matrix  of  aggregate  eigenvalues,  and  the  unitary  matrix 

Uh(-R)  =  ($(i;a  ;...;n:  R)  |^(i;  2; .  ;  n  :  R))  (18) 

provides  a  representation  of  the  Schrodinger  eigenstates  spanned  by  the  basis  of  Eq.  (15) 
in  the  form 

^(1;  2; . .  .  ;n  :  R)  =  #(1;  2; . .  .;n  :  R)  Uh(H).  (19) 

The  semi-colons  in  the  Schrodinger  eigenstates  4/  (1;  2; . . . ;  n  :  R)  are  employed  to  empha¬ 
size  that  these  solutions  are  obtained  from  Eqs.  (15)  to  (17),  in  which  the  orthonormal 
aggregate  product-function  basis  states  appearing  in  the  expression  of  Eq.  (19)  for  the  wave 
functions  are  not  individually  antisymmetric  under  all  electron-coordinate  permutations. 

4.1  Form  of  the  Hamiltonian  Matrix 

The  matrix  representative  of  the  aggregate  Hamiltonian  operator  appearing  in  Eq.  (17) 
takes  the  form45 

H (R)  =  (3>(i;  2; . . . ;  n  :  R)  |R(i,  2:5,.  , .,  n  :  R)  | 4>(i;  2; . . . ;  n  :  R)) 
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(20) 


=  V  {  H'"l  -  J2 

a=l  p=a+l 

where  H^)  and  'V<'a%!jHRnfj)  are  Hermitian  matrix  representatives  of  the  operators  of  Eq. 
(3)  and  (4)  in  the  spectral-product  basis.  The  atomic  energy  matrix  H(a)  is  a  constant 
diagonal  matrix  independent  of  atomic  position  vectors,  whereas  the  interaction  matrix 
V(a,d)  (Rap)  is  non-diagonal  and  depends  explicitly  upon  the  vector  separation  Rap  of 
the  two  indicated  atoms,  but  not  upon  the  individual  laboratory-frame  positions  Ra  and 
R/3  of  the  two  atoms,  nor  upon  the  position  vectors  of  the  other  atoms  in  the  aggregate. 
The  particularly  simple  form  of  Eq.  (20)  is  largely  a  consequence  of  the  orthogonality  of 
the  spectral-product  basis,  the  use  of  atomic  eigenstates  in  the  representation,  and  the 
atomic  pairwise-additive  nature  of  the  interaction  terms  in  the  Hamiltonian  operator  of 
Eqs.  (2)  to  (4). 

The  specific  forms  of  the  atomic-energy  and  pairwise-atomic  interaction-energy  matri¬ 
ces  of  Eq.  (20)  depend  upon  the  choice  of  ordering  rule  adopted  in  Eq.  (15).  Employing, 
for  example,  the  common  “odometer”  ordering  convention  of  letting  later  indices  run  to 
completion  prior  to  earlier  ones,  the  last  atomic  matrix  H^v)  and  the  pairwise  interaction- 
energy  matrix  \^n~1,n\Rn-in)  for  the  last  two  atoms  are  seen  to  take  block  diagonal 
forms,  with  matrices  constructed  in  the  appropriate  atomic  and  atomic-pair  product  states, 
respectively,  appearing  repeatedly  down  their  diagonals.  The  atomic-energy  and  pairwise- 
atomic  interaction-energy  matrices  of  Eq.  (20)  for  an  arbitrary  pair  of  atoms  (a,  (3)  can  be 
obtained  in  this  convention  by  starting  from  the  block  diagonal  versions  of  these  matrices 
and  performing  appropriate  row  and  column  interchanges  which  are  determined  by  the  po¬ 
sitions  of  a  and  (3  relative  to  the  last  two  indices  N  —  1  and  N,  respectively,  enumerating 
the  atoms. 

Regardless  of  the  choice  of  ordering  convention  implied  by  the  label  O  of  Eq.  (15),  the 
atomic  matrices  can  always  be  written  as  the  ordered  outer  products 

Hfo)  =  jlW  (g)  if)  <g) . . .  Effo  (g)  •  •  •  If)  }o  ,  (21) 

where 

Ef)  =  <*£*)(<  :  Ra)\H{a\i  :  Ra)\^{pa\i  :  Ra))  (22) 

is  the  diagonal  matrix  of  physical  eigenenergies  of  the  atom  a  [Eqs.  (9)  to  (11)],  arranged  by 
convention  in  energy-increasing  order  down  the  diagonal.  The  interaction-energy  matrices 
are  correspondingly  written 


(Raf>)  =  {iw  ®  i<2>  ®  ■  ■  •  v<“'»  («„/,)  0  ■  ■  ■  i;y  }o , 


(23) 


where 


v(a’d) ^Rap)  =  \i-j  :  Rap) \vW\v,j  :  Rap) |*(a-«(t;  j  :  Rap))  (24) 

is  the  interaction-energy  matrix  in  the  ordered  atomic-product  basis  for  the  (a,  (3)  pair, 
given  by  the  expression 


:  Rap)  =  {*£*)(*  :  Ra)  ®  : 


R 


o 


(25) 
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In  Eqs.  (21)  to  (25),  Ip0"1  is  the  unit  matrix  expressing  the  orthonormality  of  the 
antisymmetric  atomic  spectral  basis  [Eq.  (9)],  an  outer-matrix-product  convention  is 
adopted  for  the  multiplication  of  matrices  enclosed  in  the  brackets  {•••},  and  the  sub¬ 
script  “O”  implies  ordering  of  the  terms  in  the  energy  matrices  H(a),  'V^a’^\Rap)  and 
vectors  <I>(a’d)(j;  j  :  Rap)  in  accordance  with  the  positions  of  the  individual  functions  in 
the  row  vector  of  Eq.  (15).  The  pair-product  states  of  Eq.  (25)  are  written  for  notational 
simplicity  as  functions  of  the  vector  separation  i?ajg,  rather  than  of  the  two  position  vec¬ 
tors  Ra  and  Rg ,  in  accordance  with  primary  interest  here  in  their  use  in  evaluation  of 
V^HRaf})- 

As  should  be  clear  from  the  discussion  of  Section  3,  the  functions  3 >(a>d)(j;j  ;  Ra^ 
span  the  totally  antisymmetric  representation  of  Sna+np  once  and  only  once,  but  also  span 
other  irreducible  representations  of  the  a,  (3  atomic-pair  permutation  group,  and  so  these 
functions  do  not  warrant  a  subscript  “p”  designating  a  physical  subspace  of  the  generally 
reducible  product  of  Eq.  (25).  It  should  also  be  clear  that  the  matrix  Ep^  of  Eq.  (22)  in 
the  atomic  spectral  basis  a  has  a  dimension  smaller  than  that  of  the  matrix  v(a,P\Rap) 
of  Eq.  (24)  in  the  a1  (3  atomic-product  spectral  basis,  although  the  matrices  and 

V’(a’'8)(i?a/g)  of  Eqs.  (21)  and  (23),  respectively,  have  the  same  dimension  as  the  entire 
spectral-product  basis  of  Eq.  (15). 

The  vector  dependence  of  v(a'P\Rap)  upon  Rap  indicated  in  Eq.  (24)  arises  from 
the  dependence  of  the  product  basis  of  Eq.  (25)  on  both  the  magnitude  and  the  angular 
orientation  of  the  vector  Rap  in  the  laboratory  frame.  Using  the  transformation  properties 
of  the  atomic  basis  functions  under  rotations,65  the  interaction-energy  matrix  of  Eq.  (24) 
is  given  by  the  expression54 

^a’m(Rap)  =  (26) 

where  R(a,^)(.Ra/g)  =  {D(a)(i?a/g)  eg)  D ^  (Rap)}o  is  an  ordered  outer  product  of  Wigner 
rotation  matrices  for  each  atom  effecting  the  transformation  from  a  reference  orientation  in 
which  the  two  z  axes  of  the  interacting  atoms  are  co-aligned  to  one  in  which  they  take  their 
physical  orientations  in  the  laboratory  frame.65  The  irreducible  interaction-energy  matrix 
v(a’d)  (jRajg)  of  Eq.  (26),  which  is  obtained  from  Eq.  (24)  evaluated  in  the  co-aligned 
reference  frame,  depends  only  upon  the  scalar  separation  of  the  two  atoms  Rap ,  whereas 
the  notation  Rap  (=  <f>ap,  6ap ,  0)  employed  in  the  Wigner  functions  in  Eq.  (26)  refers  to 
the  two  angles  specifying  the  direction  of  the  vector  Rap  in  the  laboratory  frame. 

Once  the  matrices  of  Eqs.  (22)  and  (24)  have  been  constructed,  those  of  Eqs.  (21) 
and  (23)  must  be  assembled  by  performing  the  row  and  column  re-orderings  implied  by 
the  symbol  “O” .  Although  the  dimension  of  the  resulting  Hamiltonian  matrix  of  Eq.  (20) 
grows  exponentially  with  the  number  of  atoms  considered  [~  O (BN)  for  large  numbers  N 
of  identical  atoms  with  B  basis  states  each],  the  pairwise  nature  of  the  interaction  terms 
of  Eq.  (4)  insures  that  nonzero  matrix  elements  are  obtained  only  between  aggregate 
spectral-product  states  that  differ  by  no  more  than  two  atomic  eigenstates,  not  unlike  the 
familiar  Slater’s  rules  for  two-electron  atomic  matrix  elements  between  determinantal  wave 
functions  comprised  of  orthonormal  orbitals.63,66  Accordingly,  H (R)  is  a  sparse  matrix 
which  contains  non-zero  elements  the  number  of  which  in  any  row  or  column  grows  only 
as  the  square  of  the  number  of  atoms  N  and  of  the  number  of  basis  functions  B  [~ 
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0(N2 B2')\.  rather  than  exponentially,  with  large  numbers  of  atoms  and  basis  functions. 
The  latter  observation  follows  from  the  number  of  spectral-product  functions  that  can 
differ  from  a  given  spectral  product  by  zero,  one,  or  two  individual  atomic  labels,  and 
insures  that  H(i?)  grows  more  sparse  as  N  and  B  increase.  These  circumstances,  and  the 
fact  that  the  locations  of  the  non-zero  elements  in  any  row  or  column  of  the  matrix  of  Eq. 
(23)  are  known  a  priori ,  suggest  that  techniques  familiar  from  large-scale  configuration- 
interaction  methods,63  including  particularly  recursive  or  unitary  group  methods  which 
avoid  construction  of  the  entire  Hamiltonian  matrix  at  any  one  time,  can  be  employed  in 
obtaining  the  eigensolutions  of  Eq.  (17),  issues  addressed  separately  below. 


4.2  Evaluation  of  Atomic-Energy  and  Interaction-Energy  Matrices 

The  atomic-energy  matrices  of  Eq.  (22)  can  be  evaluated  once  and  for  all  employing 
conventional  electronic  structure  theory  in  its  various  forms59-63  in  constructing  the  large 
denumerable  spectrum  of  eigenfunctions  generally  required  to  satisfy  the  completeness 
conditions  of  Eq.  (16).  Of  course,  highly  accurate  experimental  values  are  also  generally 
available  for  a  significant  number  of  the  low-lying  atomic  energies  of  most  atoms.67  Accord¬ 
ingly,  adopting  a  commonly  employed  stratagem,  the  available  experimental  values,  possi¬ 
bly  corrected  to  Coulombic  limits  by  approximate  removal  of  relativistic  contributions,68,69 
can  be  employed  in  place  of  the  calculated  low-lying  values,  which  are  generally  less  ac¬ 
curate  than  the  former,  whereas  the  calculated  higher-lying  values  can  be  employed  as  a 
pseudospectrum  of  energies  in  place  of  the  actual  Rydberg  and  continuum  atomic  states.70 
Such  a  calibration  procedure,46  which  is  easily  incorporated  into  the  Hamiltonian  of  Eqs. 
(20)  to  (26),  is  particularly  required  in  cases  where  the  calculated  low-lying  atomic  ener¬ 
gies  are  not  in  accord  with  the  energy  orderings  of  the  experimental  spectra.  Employing 
this  atomic  calibration,  the  potential  energy  surfaces  obtained  from  the  development  are 
guaranteed  to  dissociate  into  accurate  atomic  separation  limits,  with  calculations  required 
only  to  provide  (Coulombic)  interaction  energies  among  the  atomic  constituents,  rather 
than  the  absolute  values  of  the  total  aggregate  energies.  Of  course,  in  many  cases  such  an 
atomic  calibration  procedure  will  be  unnecessary,  and  the  calculated  atomic  energies  can 
be  employed  without  modification  directly  in  construction  of  the  Hamiltonian  matrix. 

The  intergrals  required  in  the  irreducible  pairwise  interaction-energy  matrices  of  Eq. 
(26)  are  seen  to  be  two-center,  one-  and  two-electron  Coulombic  integrals,  which  can  all 
be  calculated  in  closed  analytical  forms  when  Gaussian39-44  or  Slater- like71, 72  orbitals  are 
employed,  ensuring  that  high  precision  can  be  achieved  in  their  evaluation.  Specifically, 
the  irreducible  interaction-energy  matrix  appearing  in  Eq.  (26)  is  written  in  the  form 


v{a^\Rap)  =  f  dra  f  dr/?  {7(a)(ra)  <8)7(/3)(r/3)}E(a,/3)(ra;r/3  :  i?a/3)  (27) 

J r^,  J  r  a 


using  Eq.  (24),  where  [cf.,  Eq.  (4)] 


V{a’P\ra-:rp  :  Raj3) 


ZaZpe2  Zpe1 


Zae2 


R, 


a  (3 


+ 


r  a  R/3 1  \rp-Ra\  \ra-rp 


(28) 
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is  the  two-center  Coulombic  interaction  operator  in  the  co-aligned  coordinate  system,  Ra  = 
Rak  and  Rp  =  Rpk,  with  k  the  unit  vector  along  the  co-aligned  z  axes,  and 


7W(r„)  =  (#<“>(*  :  Ra)\ £i(r„  -  rf )!*<“>(<  :  Ra)) 


(29) 


is  the  so-called  one-electron  transition  density  matrix  for  atom  a, 73-75  with  a  similar 
expression  employed  for  atom  /3. 

Evidently,  the  Hermitian  matrices  of  Eq.  (29),  properties  of  which  have  been  discussed 
previously  in  the  context  of  nuclear  spectroscopy  and  nucleon  scattering  by  nuclei,73-75 
determine  the  associated  irreducible  interaction-energy  matrices  of  Eqs.  (26)  and  (27), 
and  so  constitute  the  computational  invariants  required  in  the  development  of  Eqs.  (20) 
to  (26).  The  matrices  of  Eq.  (29)  formally  involve  all  the  eigenstates  of  an  individual 
atom,  and  so  they  differ  from  so-called  “reduced”  one-electron  density  matrices,  which 
are  scalar  functions  appropriate  for  individual  atomic  eigenstates.62  Rather,  the  diagonal 
elements  of  the  matrix  7^a'l(ra)  are  the  one-electron  charge-density  distributions  of  the 
spectral  eigenstates  of  atom  a,  whereas  the  off-diagonal  elements  give  the  one-electron 
transition  densities  connecting  these  eigenstates.  These  quantities  can  be  stored  in  the 
form  of  matrices  r(a)  in  the  combined  atomic-state  orbital-basis  (Ta,7a)  representation, 
and  employed  in  Eq.  (27)  in  repeated  “on-the-fly”  calculations  of  the  interaction  energies 
required  in  construction  of  the  Hamiltonian  matrix  for  the  atomic  configurations  R  consid¬ 
ered.  Alternatively,  the  pairwise  interaction-energy  matrices  of  Eq.  (27)  can  be  tabulated 
on  a  grid  of  interatomic  separations  Rap  for  the  interacting  atomic  pairs  of  interest  and 
retained  for  repeated  Hamiltonian  assembly.  Accordingly,  in  this  approach,  large-scale 
configuration-interaction  calculations  of  the  atomic  spectral  states  need  be  performed  only 
once,  and  the  asssociated  one-electron  transition-density  or  interaction-energy  matrices 
retained  for  repeated  use,  avoiding  explicit  high-level  molecular  configuration-interaction 
or  related  calculations  involving  antisymmetric  aggregate  many-electron  basis  states  at 
different  aggregate  configurations. 

4.3  Eigenspectrum  of  the  Spectral-Product  Hamiltonian  Matrix 

Although  the  development  of  Eqs.  (1)  to  (29)  can  be  employed  directly  in  computa¬ 
tional  applications  in  selected  cases,  a  number  of  factors  mitigate  against  adopting  such  an 
approach  more  generally.  Specifically,  the  spectrum  [Eq.  (1)]  of  the  Hamiltonian  operator 
of  Eqs.  (2)  to  (4)  contains  the  aforementioned  unphysical  states,  including  continua  in 
which  the  physical  solutions  are  embedded  when  any  one  atom  in  the  aggregate  contains 
three  or  more  electrons.32-37  Accordingly,  approximations  to  the  physical  states  of  Eqs. 
(1)  to  (4)  in  the  full  spectral-product  basis  may  generally  not  be  easily  discernible  among 
the  unphysical  quasi-continuum  states  also  spanned  by  the  basis,  and  may,  in  fact,  be 
computationally  inaccessible  in  the  presence  of  the  afore-mentioned  low-energy  non-Pauli 
states  present  which  are  associated  with  over-filled  inner-shell  configurations  and  their  cor¬ 
responding  continua.  This  situation  is  complicated  by  the  fact  that  precise  antisymmetry 
will  only  be  obtained  from  the  development  of  Eqs.  (1)  to  (29)  in  the  limit  of  spectral 
closure,  making  identification  of  the  desired  physical  states  potentially  subjective. 
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The  task  of  classifying  the  states  obtained  from  the  Hamiltonian  matrix  of  Eqs.  (20) 
to  (29)  is  further  complicated  by  their  degeneracies.  Specifically,  the  pairwise  interaction- 
energy  matrices  of  Eq.  (27)  will  be  diagonal  in  the  atomic  spin  quantum  numbers  Sa,  Msa 
and  Sp1Msf31  since  the  Coulomb  interaction  operator  of  Eq.  (28)  is  independent  of  spin 
coordinates.  Consequently,  the  aggregate  Hamiltonian  matrix  of  Eq.  (20)  will  generally 
factor  into  noninteracting  blocks  associated  with  particular  spin  configurations  of  the  ag¬ 
gregate  atoms  (Sa,  Msal  a  —  1,  2, ...  IV),  and  the  eigensolutions  obtained  from  Eqs.  (17)  to 
(29)  will  generally  be  fragmented  into  degenerate  orthogonal  components  associated  with 
these  contributing  spin  configurations.  These  components  must  be  appropriately  combined 
ex  post  facto  to  form  states  of  the  correct  irreducible  symmetries. 

In  a  case  in  which  only  two  different  aggregate  spin  configurations  are  present  in  the 
spectral  product  basis  [<&(i;  2 ,;n  :  R)  —  <&i(i;  2; . . . ;  n  :  R)  ©  #2(1;  2; . . . ;  N  :  it)], 
for  example,  Eq.  (17)  factors  into  the  two  independent  equations 

Hn  (it)  •  Ui  (it)  =  Ui  (it)  •  Ei  (it)  (30a) 

H 22(H)  •  U2(it)  =  U2(it)  •  E2(it),  (306) 

where  the  off-diagonal  matrix  elements  H12(it)  and  H21(it)  vanish,  and  the  diagonal  en¬ 
ergy  matrices  Ei  (it),  E2(it)  and  solution  matrices  Ui  (it),  U2(it)  are  obtained  from  the 
indicated  separate  diagonalizations  of  Hu  (it)  and  H22(it),  respectively.  The  eigenfunc¬ 
tions  of  Eq.  (19)  in  this  case, 


^1(1;  2; . 

•  ;n  :  R)  =  $1(1;  2;. 

.  .;n:  R)  -Ui(Jl) 

(31a) 

■  :n  :  R)  =  #2(i;2;.. 

.\n:  R)  •  U2(il) j 

(316) 

may  include  a  solution  [M>i(i;  2; . . . ;  n  :  R)]k  which  in  the  limit  of  closure  has  the  same 
eigenvalue  as  a  solution  [\l>2(i;  2; . . . ;  n  :  R)]i  if  both  spin  configurations  can  contribute  to 
the  total  wave  function  of  the  chemical  aggregate  [{Ei (it)}fc&  ~  {E2(it)};;].  In  general, 
these  two  components  must  be  identified  and  appropriately  combined  to  construct  the 
correct  antisymmetric  aggregate  eigenfunction  which  together  they  comprise.  The  two  de¬ 
generate  spin  multiplets  obtained  from  Eqs.  (31),  however,  need  not  necessarily  transform 
irreducibly  under  electron  coordinate-label  permutations,  and  so  there  may  be  no  way  in 
which  to  readily  identify  them  a  priori  as  potentially  degenerate  components  of  a  totally 
antisymmetric  physical  eigensolution  of  the  Schrodinger  equation,  particularly  in  the  pres¬ 
ence  of  other  nearly  degenerate  unphysical  states.  Although  these  circumstances  may  seem 
unfamiliar,  a  similar  situation  would  arise  in  treatments  of  atomic  structure  employing 
unrestricted  spin-orbital  product  configurations  in  constructing  the  Hamiltonian  matrix, 
rather  than  Slater  determinants.66  Combinations  of  products  of  spatial  functions  and  spin 
eigenstates,76,77  or  functions  of  spatial  coordinates  only,78  also  employed  in  constructing 
eigensolutions  of  irreducible  symmetry  similarly  correspond  to  solutions  expressed  in  the 
fragmented  forms  of  Eqs.  (30)  and  (31). 

Finally,  the  factoring  of  the  spectral-product  Hamiltonian  matrix  into  individual  atomic 
spin  blocks  can  also  lead  to  a  failure  to  provide  the  correct  numbers  of  multiplet  states  in 
limits  of  large  interatomic  separations.  This  will  be  so  in  the  common  case  of  dissociations 
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into  open-shell  atoms,  where  the  failure  to  achieve  closure  in  the  spectral-product  basis 
in  the  absence  of  explicit  electron  exchange  terms  can  give  results  which  include  only  the 
lowest-lying  spin  multiplet  of  an  allowable  number  of  degenerate  states  at  large  atomic 
separations.  These  circumstances  can  lead  to  calculations  of  highly  unphysical  potential 
energy  surfaces  which  correlate  with  incorrect  dissociation  limits.  Similar  remarks  apply 
to  cases  in  which  charge-transfer  terms  contribute  significantly  to  eigensolutions  at  large 
interatomic  separations,  where  the  failure  to  achieve  closure  in  the  spectral-product  ba¬ 
sis  and  in  the  corresponding  development  of  Eqs.  (17)  to  (29)  could  introduce  spurious 
variations  in  potential  energy  surfaces  with  increasing  interatomic  separations. 

5.0  Convergence  in  the  Spectral-Product  Basis 

Group-theoretical  methods  have  been  used  previously  in  understanding  the  conse¬ 
quences  of  the  change  in  electron-coordinate  permutation  symmetry  group  upon  formation 
of  atomic  aggregates,32  and  have  been  employed  above  to  verify  the  presence  of  the  totally 
antisymmetric  representation  in  the  spectral-product  basis.  Such  methods,  however,  are 
largely  inapplicable  to  reduction  of  the  basis  in  the  group  Snt  for  computational  purposes 
in  light  of  the  non-symmetrical  nature  of  the  representation  and  the  associated  absence 
of  all  electron  coordinate  permutations  in  explicit  form  in  Eq.  (15).  In  the  present  Sec¬ 
tion,  the  totally  antisymmetic  representation  of  Snt  spanned  by  the  basis  of  Eq.  (15)  is 
isolated,  the  additional  troublesome  issues  identified  in  the  immediately  preceding  Section 
are  overcome,  and  convergence  to  the  correct  Schrodinger  eigenstates  in  the  antisymmetric 
subspace  of  the  spectral-product  representation  is  demonstrated. 

Symmetry-adapted  eigensolutions  are  devised  in  the  following  in  the  blocked  forms 


*(i;2;...;IV:R)  =  {^)(i;  2; . . . ;  M  :  H),  ^(i;  a; . . . ;  IV  :  H)}, 


(32a) 


E  (R) 


( V{i\R)  0  \ 
V  0  E  (su)(R)J 


(32  b) 


where  ^^(1;  2; . . . ;  n  :  R)  and  \l/(w)(i;  2; . . . ;  n  :  R)  refer  to  physical  “p”  and  unphysical 
“it”  solutions  of  Eqs.  (17)  to  (19),  respectively,  and  Eg^(R),  EgU)(R)  are  the  corresponding 
diagonal  eigenvalue  matrices.  These  eigenspectra  are  obtained  formally  from  a  unitary 
transformation  of  the  aggregate  Hamiltonian  matrix  of  Eq.  (20), 


Hs(fl)  =  U8(J1)'  ■  Hlfij  ■  Uslft)  -> 


(H  S\R) 
V  0 


(33) 


which  explicitly  isolates  the  physical  Hg^  (R)  and  unphysical  HgU')  (R)  blocks  of  the  spectral- 
product  Hamiltonian  matrix.  The  unitary  transformation  matrix  Us(R)  obtained  from 
diagonalization  of  the  aggregate  metric  matrix  S  (R)  constructed  from  antisymmetrized 
spectral-product  states  is  shown  to  accomplish  the  blocking  of  Eqs.  (33),  and  the  corre¬ 
sponding  eigenfunctions  and  eigenvalues  of  Eqs.  (32)  are  shown  to  converge  in  the  limit  of 
closure  to  results  obtained  from  standard  variational  developments  employing  conventional 
prior  antisymmetry  constraints. 

5.1  Prior  Antisymmetry 
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The  development  leading  to  Eqs.  (32)  and  (33)  begins  with  construction  of  an  antisym¬ 
metric  spectral-product  basis,  familiar  from  the  early  development  of  Moffitt,46  obtained 
employing  the  prior  projection 

#(i,  2, . . . ,  n  :  R)  =  PA(i,  a, . . . ,  n)#(i;  2; . . . ;  n  :  R),  (34) 

where  d?(i;  2; . . . ;  n  :  R)  can  be  taken  to  be  a  subspace  of  Eq.  (15)  having  good  total  Ms 
compatible  with  the  total  S  value  of  interest  (Ms  <  S ).  The  projector 

nt\ 

PA(  1,  2, . .  ,,n)  =  (m \n2\  ■  ■■nN\)~1/2(ntl)~1/2  J^(-l )6pPp  (35) 

p=i 

is  the  familiar  nt-electron  antisymmetrizer  “wave-function  normalized”  in  the  limit  R  — > 
00. 63  In  Eq.  (35),  5p  is  the  parity  of  the  permutation  Pp,  the  sum  is  over  all  nt\  electron 
coordinate  permutations,  and  the  correction  factor  (nd^!  •  •  -n.jv!)-1//2  insures  asymptotic 
wave-function  normalization  in  the  presence  of  the  prior  antisymmetry  of  the  individual 
atomic  states  in  Eq.  (15).  The  use  of  commas  between  electron  coordinates  on  the  left- 
hand  side  of  Eq.  (34)  serves  to  indicate  that  the  individual  terms  in  the  basis  are  of  good 
permutation  symmetry  type,  as  opposed  to  those  of  the  spectral  product  basis  of  Eq.  (15). 

Employing  Eqs.  (15)  and  (16),  the  basis  of  Eq.  (34)  is  seen  to  satisfy  the  closure 
relation 


<I?(i,  2, . . . ,  n  :  R)  •  $(i/,  2 n'  :  R )^ 


ntl 


N 


n\\n2\ 


■nN\ 


n  SfHi-i')  (36) 


a= 1 


in  the  totally  antisymmetric  subspace  of  Eq.  (15).  Consequently,  the  antisymmetrized 
spectral-product  or  Moffitt  basis  is  formally  Q-fold  redundant  [Q  =  nt\/(ni\ri2].  ■  ■  •  njvO] 
for  representations  of  antisymmetric  states.25,27-29  To  clarify  the  origins  of  this  circum¬ 
stance,  note  that  although  PA  of  Eq.  (35)  nulls  the  unphysical  irreducible  representations 
of  Snt  contained  in  the  reducible  spectral-product  basis,  the  failure  to  exclude  the  re¬ 
dundant  contributions  that  arise  in  antisymmetric  states  formed  from  unrestricted  orbital 
products  in  the  spectral- product  basis  gives  rise  to  the  factor  Q  in  Eq.  (36).  Such  over¬ 
counting  is  avoided  within  each  individual  set  of  atomic  states  (i  :  Ra )  by  including  all 
orderings  of  spin-orbital  configurations  but  excluding  explicit  permutations  among  electron 
coordinates  i,  or,  equivalently,  by  limiting  the  contributing  spin-orbital  products  to  strictly 
ordered  electron  configurations  in  Slater  determinants.62  However,  since  no  restrictions  are 
placed  on  the  contributions  in  Eq.  (34)  arising  from  products  of  spin-orbitals  from  two 
or  more  different  atomic  basis  sets,  and  all  electron  coordinate  permutations  among  these 
are  included  implicitly  in  the  limit  of  closure,  Eq.  (34)  contains  redundant  terms  in  this 
limit.25,27-29,73-75  More  specifically,  the  factor  nf.  in  the  numerator  of  the  redundancy 
factor  Q  is  the  total  number  of  possible  arrangements  of  nt  electrons  among  nt  available 
positions,  whereas  the  factor  •  •  •  un]-  in  the  denominator  of  Q  corrects  this  value 

only  for  re-arrangements  within  the  individual  atoms,  which  do  not  arise  in  forming  the 
linearly  independent  atomic  spectral  states  by  using  a  single  assignment  of  electrons  to 
orbitals  or  the  equivalent  restriction  to  ordered  configurations  in  Slater  determinants.62 
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The  overcompleteness  of  Eq.  (36)  is  of  no  consequence  in  the  limit  R  — >  oo,  since  the 
interaction  terms  in  the  Hamiltonian  matrix  vanish  identically  and  the  Hamiltonian  matrix 
correspondingly  has  no  off-diagonal  elements.  In  this  case,  the  states  of  Eqs.  (15)  and  (34) 
give  the  same  aggregate  energies,  and  the  latter  otherwise  provides  results  equivalent  to  use 
of  the  former  atomic-product  states.  Moreover,  for  large  values  of  R.  when  perturbation 
theory  can  be  employed,  the  factor  Q  can  be  taken  into  account  explicitly  in  expressions 
for  wave  functions  and  energies  involving  spectral  summations,  although  care  must  be 
exercised  in  this  connection.25-38  More  generally,  however,  for  arbitrary  finite  values  of 
-R,  where  the  basis  of  Eq.  (34)  is  no  longer  orthonormal  and  perturbation  theory  is 
unsuitable,  the  redundancy  in  the  limit  of  closure  gives  rise  to  linear  dependence  which 
must  be  overcome  in  computational  applications.  As  will  be  demonstrated,  it  is  possible 
to  not  only  employ  familiar  canonical  orthogonalization  procedures  to  remove  the  linear 
dependence  in  Eq.  (34), 79  but,  correspondingly,  to  also  remove  the  unwanted  non-Pauli 
representations  of  the  symmetric  group  Snt  present  in  the  spectral-product  basis  of  Eq. 
(15). 

5.2  Metrically-Defined  Hamiltonian  Representation 

Following  standard  procedures,79  the  potentially  overcomplete  Mofhtt  basis  of  Eq.  (34) 
is  employed  in  a  linear  variational  solution  of  the  aggregate  Schrodinger  equation,  giving 
the  familiar  matrix  equation 

H(R)  ■  Uh{R)  —  S(R)  ■  Uh(R)  •  E(R)  (37) 

appropriate  in  a  non-orthonormal  representation.63  Here, 

N  N 

H(R)  =  {  iT(a)(R)  +  y(a,/3)(#)}  (38) 

a= 1  p=a+ 1 

is  the  metrically-defined  Hamiltonian  matrix  in  the  totally  antisymmetric  basis,  where  the 
individual  atomic  H^a\R)  and  atomic-interaction  V^a,^(R)  terms,  corresponding  to  the 
operators  of  Eqs.  (3)  and  (4),  respectively,  are  now  dependent  upon  the  positions  of  all 
the  atoms  in  the  aggregate  consequent  of  this  antisymmetry  [cf.,  Eq.  (20)], 

S (R)  —  (<l?(i,  2, . . . ,  n  :  R) |3>(i,  2, . . . ,  n  :  R))  (39) 

is  the  overlap  or  metric  matrix  of  the  antisymmetric  basis,  E(R)  is  the  aggregate  diagonal 
energy  eigenvalue  matrix,  and  Uh{R )  is  the  solution  matrix  that  contains  the  column 
eigenvector  representations  of  the  antisymmetric  Schrodinger  states  in  the  non-orthogonal 
basis  of  Eq.  (34), 


^(i,  2, . . . ,  n  :  R)  =  $(i,  2, . . . ,  n  :  R)  ■  UH(R)-  (40) 

The  matrices  S  (R)  and  Uh{R )  together  determine  the  projection  of  the  basis  states 
of  Eq.  (34)  on  the  solutions  of  Eq.  (40)  in  the  usual  form, 

(#(i,  2, . . . ,  n  :  R)  |*(i,  2, . . . ,  n  :  R))  =  S(R)  ■  UH(R),  (41) 
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where  Uh(R )  satisfies  the  familiar  modified  orthogonality  condition 

Uh(RY  ■  S(R)  -Uh(R)  =  I,  (42) 

and  I  is  the  unit  matrix  of  dimension  equal  to  that  of  the  basis  of  Eq.  (34).  It  should  be 
noted  that  different  fonts  are  employed  in  Eqs.  (17)  to  (20)  and  in  Eqs.  (37)  to  (42)  to 
distinguish  between  matrices  and  other  quantities  represented  in  the  spectral-product  and 
Moffitt  basis  sets,  respectively. 

5.3  Removal  of  Linear  Dependence  in  the  Moffitt  Basis 

In  the  closure  limit  of  Eq.  (36),  the  Moffitt  basis  of  Eq.  (34)  is  linearly  dependent,  the 
metric  matrix  of  Eq.  (39)  correspondingly  becomes  singular  [det|S(i?)|  —¥  0],  and  Eq.  (37) 
is  generally  ill  conditioned.  The  linear  dependence  of  the  basis  of  Eq.  (34)  in  the  limit 
of  closure  is  avoided  by  isolating  a  linearly  independent  subspace  employing  the  unitary 
transformation  matrix  that  diagonalizes  S  (i?), 

Us(R)t  •  S (R)  ■  Us (R)  =  s d(R),  (43) 

where  the  diagonal  matrix  s d{R)  contains  a  positive  block  [s^(JR)]  and  a  block  [s^(il)] 
that  tends  to  zero  value  in  the  closure  limit.  By  appropriate  orderings  of  the  eigenvalues  and 
vectors  of  S(R),  the  non-zero  block  associated  with  the  column  eigenvectors  in  Us(-R)  that 
provide  the  physically  significant  linearly-independent  subspace  of  Eq.  (34)  can  be  placed 
in  the  upper  left-hand  corner  of  s d(R),  whereas  the  block  that  tends  to  zero,  associated 
with  the  linearly-dependent  or  “unwanted”  portion  of  the  basis  of  Eq.  (34) ,  can  be  placed 
in  the  lower  right-hand  corner,  giving 


s  d(R) 


(^’(R)  o  V  ./Qip  °1 
1  o  S<”W  Vo  o  J- 


(44) 


Here,  Ip  is  the  unit  matrix  in  the  linearly-independent  physical  subspace,  0(t  is  the  zero 
matrix  in  the  unwanted  subspace,  and  Q  is  the  redundancy  factor  of  Eq.  (36).  The  right- 
hand  side  of  Eq.  (44)  follows  from  the  observation  that  the  eigenvalues  of  the  metric 
matrix  S(R)  are  either  Q  or  0  in  the  limit  of  closure,  in  accordance  with  the  presence  of 
this  factor  in  Eq.  (36)  and  the  linear  dependence  of  the  unwanted  subspace. 

Employing  the  matrix  Us(-R)  of  Eq.  (43),  the  basis  of  Eq.  (34)  is  transformed  into  a 
new  basis  $s(15  2, . . . ,  n  :  R)  comprising  the  linearly-independent  subspace  $g^(i,  2, . . . ,  n 
R)  in  which  the  physical  eigenstates  can  be  constructed,  and  the  unwanted  subspace 
<l?g^(i,  2, . . .  ,n  :  R).  The  two  subspaces  are  given  by  the  common  expression 

$s(i,  2,  ,  n  :  R)  =  $(i,  2,  ,  n  :  R)  ■  VS(R)  ■  s d(R)~1/2 

=  a, . . . ,  n  :  R),  ^^(u  2, . . . ,  n  :  R)},  (45) 

where 

*£°(i,  2,  -  -  - ,  n  :  R)  =  {$(i,  2, . ,  n  :  R)  ■  US(R)}P  •  sJj)  (R)-1/2  (46) 
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is  the  orthonormal  linear ly-independent  subspace  and  {. .  .}p  designates  the  physical  part 
of  the  enclosed  row  vector.  The  linearly-dependent  subspace 

*£°(i,  a,  •  •  • ,  n  :  R)  =  {*(1,  a, . . . ,  n  :  R)  •  Us (R)}u  •  sj0  (i*)”1/2  (47) 

is  ill-defined  in  the  limit  (R)  — >■  0r  but  is  not  required  here.  The  closure  of  the  physical 

basis  for  antisymmetric  states, 

N':R)1  = 

{$(1,  2, . . . ,  N  :  R)  •  Us  (R)}p  •  s^GR)-1  •  {*(i',  2 N' :R)-  U  s(i?)}+  -> 

N 

Q-'m  1,  a, . . . ,  JV  :  H)  •  *(i',  a',.. . ,  AT'  :  B)}„  ->  {  n  (48) 

a=l 

follows  from  Eqs.  (36)  and  (44)  and  the  unitarity  of  Us(-R),  with  the  redundancy  factors  Q 
that  appear  in  the  numerator  and  denominator  of  the  left-hand  side  of  the  last  line  in  Eq. 
(48)  canceling  in  the  limit  of  completeness.  The  second  subscript  “p”  on  the  right-hand  side 
of  Eq.  (48)  refers  to  restriction  of  the  delta- function  product  to  the  linearly-independent 
subspace. 

The  matrix  Us(-R)  of  Eq.  (43)  can  be  employed  in  the  usual  way  to  transform  the 
potentially  ill-conditioned  Hamiltonian  matrix  H(R)  of  Eqs.  (37)  and  (38)  to  a  form 
that  isolates  its  linearly-independent  block,  giving  the  generally  well-conditioned  matrix 
Schrodinger  equation 


H^(R)  ■  U%\r)  =  U%\R)  ■  E^(R),  (49) 

where  the  reduced-dimension  matrix 

H^(R)  =  3, . . . ,  N  :  fl)|H(l,  2, . . . ,  JV  :  R) 2, ....  i V  :  «)) 

=  af( R)~U2  ■  {Us(B)*  ■  H(R)  ■  Us(B)}pp  •  4ri(R)~1/2 
->  Q-1/2{Us(R)t  ■  H(R)  ■  Us(R)}ppQ-,/2 

=  E{Hs“’(«)+  E  ^“^(fi)}  (50) 

a=l  /3=a+l 

provides  the  non-singular  portion  of  the  metrically-defined  Hamiltonian  matrix.  Here,  the 
notation  {. .  .}pp  implies  that  only  the  upper- left  (pp)  block  of  the  enclosed  matrix  is  to  be 
retained,  and  the  individual  atomic 

nia)(R)  =  s(f\R)-1/2  ■  {Us(B)f  -  H(“>(B)  •  Us(R)}pP  ■  sP(B)-1'2 

-4  Q~l/2  {Us(fl)*  •  JTM(R)  •  U s(B)}pp  Q-1'2  (51) 
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and  atomic-interaction 

V^"'fj)(R)  =  s {f\R)-1/2  •  (Us(-R)^  •  V^\R)  •  Us(iZ) }pP  ■  s {f\R)~1/2 

Q~1/2  {Us(^)t  •  V^\R)  •  Us(R)}PP  Q~1/2  (52) 

terms  in  Eq.  (50)  refer  to  transformations  of  the  corresponding  terms  in  Eq.  (38). 

The  matrix  of  physical  energies  (R)  and  of  corresponding  column  eigenvectors 
ug\R)  of  Eq.  (49)  are  of  reduced  dimensions  relative  to  E(R)  and  Uh(R)  of  Eq.  (37), 

respectively.  Since  the  eigenstates  $gP')  (l,  2, . . .  ,  n  :  R)  of  the  S (R)  matrix  of  Eq.  (41)  are 
symmetry  adapted,  the  Hamiltonian  matrix  of  Eq.  (50)  will  be  block  diagonal  in  the  total 
spin  quantum  numbers  and  in  other  appropriate  aggregate  symmetry  labels,  in  accordance 
with  Schur’s  lemmas.60  Of  course,  the  limits  (R)  — *  Q  lp  and  s ^  (R)  — >  0U  of  Eq. 
(44)  apply  only  when  closure  is  achieved,  with  the  partitioning  of  s d(R)  into  s(f\R)  and 

s ^  (R)  somewhat  subjective,  and  the  diagonal  matrix  (R)  appearing  in  Eqs.  (50) 
to  (52)  required  more  generally,  rather  than  its  limiting  value  Q  Ip.  As  a  consequence, 
selection  of  the  dimension  of  the  physical  block  {. .  .}PP  of  the  transformed  Hamiltonian 
matrix,  which  will  generally  have  non-zero  off-diagonal  blocks  connecting  the  “p”  and  “u” 
subspaces  in  finite  representations,  is  correpondingly  subjective,  and  must  be  determined 
in  the  course  of  calculations  on  basis  of  the  accuracy  of  results  desired,  with  convergence 
achieved  when  the  off-diagonal  terms  are  judged  to  be  negligible. 

5.4  Isolation  of  the  Antisymmetric  Spectral-Product  Subspace 

The  foregoing  standard  development  employed  to  remove  linear  dependence  in  the 
prior  antisymmetrized  or  Mofhtt  basis  of  Eq.  (34)  can  also  be  employed  to  separate  the 
physical  and  unphysical  subspaces  of  the  spectral-product  representation,  and  to  construct 
corresponding  energies  and  eigenfunctions  of  the  Hamiltonian  matrix  of  Eq.  (17)  which 
converge  to  the  symmetry- adapted  eigenstates  obtained  from  the  development  of  Eqs.  (37) 
to  (52).  Specifically,  the  important  identity 

($(i;  2 ; . .  .;n  :  R)\PA(x,2, . .  .,n)|$(i;a; . .  .;n  :  R))  =  Q~lf 2  S (R),  (53) 


which  follows  from  substitution  of  Eq.  (34)  into  Eq.  (39)  and  is  valid  in  any  denumerable 
basis,  indicates  that  the  matrix  representative  of  the  antisymmetric  projector  of  Eq.  (35) 
in  the  spectral-product  basis  is  proportional  to  the  metric  matrix  S  (R)  constructed  in 
the  antisymmetric  basis  4>(i,  2, . . . ,  n  :  R).  Accordingly,  diagonalization  of  the  matrix 
S(R)  in  Eq.  (43)  is  equivalent  to  diagonalization  of  the  matrix  representative  of  the 
antisymmetric  projector  in  the  spectral-product  basis,  and  to  the  construction  of  a  subspace 
of  antisymmetric  eigenstates  in  this  representation. 

Applying  the  transformation  Us  {R)  to  Eq.  (53)  gives 


Us(H)f  ■  <*(i;  a; . . n  :  R)\Pa(i,  a, . . . ,  n)|S(i ;  a; fl))  •  Us(fl) 

=  Q- 1/2  f*T  («)  0 


o  s{;\r)J’ 


(54) 
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where  Eq.  (43)  has  been  employed  on  the  right-hand  side.  Defining  the  transformed  basis 
$s(i;2 R)  =  $(i;2 R)  U S(R) 

(55) 


=  {^«p)(i;2; ..  ,;n  :  R),  $£u)(i;  2; . . . ;  n  :  R)}, 


where 


#^/w)(i;2;..,;n  :  R)  =  {$(1;  2; . . . ;  n  :  R)  U  s(-R)}P/«  (56) 

correspond  to  the  physical  and  unphysical  subspaces,  it  follows  that  Eq.  (54)  can  be 
written  as  the  two  equations 


(*s°(i;  2;  •  •  • ;  n  :  R)  1^  ^(-l)^Pp|#^(1;  a; . . . ;  n  :  R))  =  s {p\R)/Q  Ip  (57 a) 

Tit' 


nt\ 


P= 1 


$Lu)(i;2;. 


1 


nt\ 


■ ; "  :  H)  I A  . . . ;  n  :  R))  =  s?’(R)/Q 

lit*  I 
P=  1 


.(«)  , 


0, 


(576) 


Equation  (57a)  demonstrates,  independent  of  group-theoretical  considerations,60  that  the 
“p”  subspace  of  Eqs.  (55)  and  (56)  contains  the  totally  antisymmetric  representation  of 
Snt  only  once  [s('f\R)/Q  — >  Ip],  whereas  Eq.  (57b)  demonstrates  that  the  “it”  subspace 
consists  of  non-totally-antisymmetric  eigenstates  in  the  limit  of  closure  [s^\R)/Q  — »  0r]. 
Note  that  the  factor  (1  /nt\)  employed  in  Eqs.  (57)  is  appropriate  for  a  normalized  Young 
operator,  rather  than  a  normalized  wave  function  [cf.,  Eq.  (35)].  Finally,  from  Eqs.  (16) 
and  (56)  and  the  unitarity  of  Us(-R), 


(i;a; ..  ,;n  :  R)  ■  &(sp/u)(i';2 n'  :  R)f  -»■  {  JJ  5^a)(i  -  i')}p/t 


N 


(58) 


a=l 


expresses  the  closure  of  both  subspaces  of  the  eigenstates  of  the  totally  antisymmetric 
projector  in  the  spectral-product  basis. 

It  is  seen  from  the  development  of  Eqs.  (53)  to  (58)  that  both  “p”  and  “it”  subspaces 
of  Eqs.  (55)  and  (56)  can  be  constructed  employing  the  transformation  matrix  of  Eq.  (43) 
obtained  from  the  metric  matrix,  with  “p”  referring  to  the  totally  antisymmetric  physical 
subspace  of  the  spectral-product  basis  and  “it”  to  a  mathematically  well-defined  orthonor¬ 
mal  unphysical  subspace  in  this  representation.  It  may  be  concluded,  therefore,  that  the 
transformation  matrix  Us(-R)  of  Eq.  (43),  defined  to  eliminate  linear  dependence  in  the 
explicitly  antisymmetrized  Moffitt  basis  <&(i,  2, . . . ,  n  :  R)  of  Eq.  (34),  correspondingly 
separates  the  totally  antisymmetric  and  non-Pauli  subspaces  of  the  spectral-product  basis 
3>(i;  2; . . . ;  n  :  R)  of  Eq.  (15).  Accordingly,  Us(-R)  can  be  employed  to  partition  the 
Hamiltonian  matrix  of  Eq.  (20)  into  its  physical  and  unphysical  blocks  [Eq.  (33)],  and  to 
construct  physical  Schrodinger  eigenstates  in  the  spectral-product  basis  [Eqs.  (32)]. 

Employing  the  transformation  matrix  Us(-R)  of  Eq.  (43),  the  physical  block  of  the 
matrix  Schrodinger  equation  in  the  spectral-product  basis  [Eq.  (17)]  takes  the  form  [cf., 
Eq.  (49)] 

■  Ug>(fl)  =  U^(fl)  ■  E<f (59) 
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where 


H {g\R)  =  2; . . . ;  N  :  R)\H(i,  2, . . . ,  N  :  R) (1;  2; . .  $N  :  R)) 


-  /*(p) 


(p). 


AT 


iV 


=  {U^(H)-H(H).Us(H)}w  =  X){Hs“)(«)+  E  Vs“’W(«)} 

a=l  /3=a+l 

is  the  physical  block  of  the  Hamiltonian  matrix  of  Eq.  (20),  with 


(60) 


H =  {US(H)'  ■  H<“>  ■  Us(fl)}w, 


(61) 


and 

V{sa’P)(R)  =  {Us («)t  •  V{a’P)(Rap)  ■  Us (R)}PP,  (62) 

corresponding  to  transformations  of  the  individual  terms  of  Eqs.  (21)  and  (23),  now 
dependent  upon  the  positions  of  all  the  atoms  in  the  aggregate  as  a  consequence  of  the 
non-local  effects  of  antisymmetry. 

5.5  Equivalence  of  Prior  and  Post  Antisymmetrization 

To  demonstrate  the  equivalence  of  results  obtained  from  the  post  and  prior  antisym¬ 
metrization  procedures  described  in  the  foregoing  developments,  note  from  Eq.  (34)  that 
the  states  <frg^(i,  2, . . . ,  n  :  R)  and  $^(1;  2; . . . ;  n  :  R)  of  Eqs.  (46)  and  (56)  are  identical 
in  the  limit  of  closure, 


4P)(U  2,  •  ■  ■ ,  n  :  R)  =  £4*^(1;  2; . . . ;  n  :  R)  •  s? ){R)-1/ 2 


,(p) 


(p), 


ntl 

s)(l;a;---;n:  R)  *s’)(l;a;...;n  :  R),  (63) 

p= 1 

where  the  left-hand  side  of  the  last  line  is  a  consequence  of  Eq.  (44)  and  the  right-hand 
side  of  the  antisymmetry  of  the  basis  <&g^(i;  2; . . .  ;n  :  R )  in  the  limit  of  closure,  the  nt\ 
permutations  of  the  unnormalized  antisymmetrizer  Ylp=i(~ ^-YpPp  providing  nt !  identical 
terms  which  cancel  the  corresponding  factor  in  the  denominator  in  this  limit.  Alternatively, 
employing  Eqs.  (53)  to  (57),  the  spectral  overlap  of  the  two  sets  of  basis  states  is  given  by 
the  expression 


(^sV;2;  ••  -  ;n  :  -R)|^s)(1,2,  -  ,™  :  R))  =  Q  1/2  (R)1/2  ^  Ip,  (64) 


in  accordance  with  these  basis  states  becoming  term-by-term  identical  functions  in  the 
closure  limit  of  Eq.  (63). 

The  foregoing  results  can  also  be  expressed  in  terms  of  the  representation  of  the  trans¬ 
formed  Mofhtt  basis  [Eq.  (46)]  in  the  transformed  spectral-product  basis  [Eq.  (56)], 


$s(i,2, . .  ,.,n  :  R)  =  $s(i;a;  ...;n  :  R)  ■  ($s(i;a; . .  ;n  :  R)  |$s(i,  2, . . . ,  AT  :  R))  -> 
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$s(i;2;...;n:  R)  Q  1/2  sd(R)  ■  sd(R)  1/2  ->  $s(i;  2;  •  •  •  5  n  :  R)  ■  ^Ip  ®  j  ,  (65) 

employing  Eqs.  (53)  to  (57)  and  the  closure  limit  in  the  last  line.  It  is  seen  that  Eq. 
(65)  corresponds  to  Eq.  (63)  for  the  physical  subspaces  of  the  two  basis  sets,  and  also 
exhibits  explicitly  the  linear  dependence  of  the  subspace  3>g^(i,  2, . . . ,  n  :  R)  in  the  limit 
of  closure. 

Similarly,  introducing  the  closure  relation  of  Eq.  (58)  for  the  physical  states  on  either 
side  of  the  Hamiltonian  operator  in  Eq.  (50),  and  using  Eq.  (64),  gives 

=  Q~1/2  s {f\R)1/2  •  H {g\R)  •  s {f\R)1/2  Q~1/2  H {g\R),  (66) 

employing  the  limiting  expression  s^\r)  — *  QIP  in  the  last  line,  where  Hg’^(Jl)  is  given 
by  Eq.  (60).  The  individual  atomic 

H{sa)(R )  =  Q-1/2  S{£\R)+1/2  ■  ■  S{£\R)+1/2  Q~1/2  ->•  H £\R)  (67) 

and  pair-interaction 

V^\R)  =  Q~1/2  S{f\R)+1/2  ■  V^(il)  •  S{f\R)+1/2  Q~1/2  ->  V^j/3)(i?),  (68) 

terms  of  Eq.  (50)  are  also  seen  to  be  equivalent  to  the  corresponding  terms  Hga')  (R)  and 

Vg*’^)  (R)  appearing  in  Eqs.  (60)  to  (62)  in  the  limit  of  closure.  Equations  (66)  to  (68) 

show  that  the  Hamiltonian  matrices  of  Eqs.  (50)  and  (60),  and  the  associated  individual 
atomic-energy  and  atomic  pairwise  interaction-energy  matrices  are  related  by  the  indicated 
renormalizations  in  finite  basis  sets,  whereas  they  are  identical  in  the  closure  limit. 

Finally,  it  follows  from  the  foregoing  development  that  the  eigenvalues  and  eigenvectors 
obtained  from  Eq.  (59)  converge  to  those  of  Eq.  (49)  [E^ (R)  — >  Eg ^(-R),  U^\R)  — > 
U^(il)],  and,  in  view  of  the  equivalence  of  the  two  basis  sets  [Eqs.  (63)  to  (65)],  that  the 
eigenstates  obtained  from  the  two  developments, 

^(p)(i,2,  ...,n  :  R)  =  ^p)(i,2,...,n  :  R)  ■  (R)  (69  a) 

*(p)(i;  2; . . . ;  n  :  R)  =  <^p)(i;  2; . . . ;  n  :  R)  •  U g}(fl),  (69 b) 

provide  identical  antisymmetric  solutions  of  the  Schrodinger  equation  in  the  limit  of  closure. 

6.0  Computational  Implementation  and  Illustrative  Calculations 

To  demonstrate  that  the  new  approach  described  in  the  foregoing  Sections  is  computa¬ 
tionally  viable,  illustrative  calculations  of  the  lowest-lying  singlet  and  triplet  states  of  the 
electron  pair  bond  (H2)  are  given  as  an  example  of  the  nature  of  the  metric  matrix  of  Eq. 
(39)  and  its  eigenvalues,  and  of  the  convergence  achieved  employing  standard  Gaussian 
basis-set  methodology.39,40  A  computationally  efficient  implementation  of  the  theory  is 
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described  which  employs  recursive  construction80,81  of  invariant  subspaces  of  the  spectral 
basis  starting  from  appropriate  reference  states  in  calculations  of  potential-energy  curves. 
Construction  of  the  entire  metric  matrix  is  avoided  in  this  way,  and  only  individual  rows 
of  the  spectral-product  Hamiltonian  matrix  need  be  evaluated  sequentially,  avoiding  its 
complete  construction  and  storage  at  any  one  time. 

6.1  Illustrative  Calculations  -  The  Electron  Pair  Bond 

In  the  two-electron  case  (nt— 2),  the  spin  functions  are  factored  out  and  attention  fo¬ 
cused  on  spatial  functions  which  are  symmetric  or  antisymmetric  under  electron  transpo¬ 
sition.  These  provide  a  basis  for  irreducible  representations  of  S2  corresponding  to  single- 
and  triplet-state  symmetries,  respectively,  with  the  singlet  states  playing  the  role  of  the 
unphysical  subspace  in  this  example.  That  is,  although  the  spatially  symmetric  states  are 
physically  significant  for  H2,  they  nevertheless  represent  an  unwanted  symmetry  type  in 
the  context  of  the  development  of  Section  5  relative  to  the  spatially  antisymmetric  states, 
and  so  can  be  used  to  illustrate  the  permutation-symmetry  isolation  method  described 
there.  In  the  absence  of  other  irreducible  representations  for  nt— 2,  both  sets  of  solutions 
can  also  be  constructed  by  direct  diagonalization  of  the  spectral-product  Hamiltonian  ma¬ 
trix  for  comparison  with  results  obtained  from  the  unitary  transformation  of  Eqs.  (59)  to 
(62)  for  isolating  the  subspace  representations. 

The  spectral-product  basis  of  Eq.  (15)  in  this  case, 

*(1;  2  :  R)  =  {<£(a)(i  :  Ra)  <g>  0^(a  :  Rp)}  q  ,  (70) 

corresponds  formally  to  all  products  of  all  discrete  and  continuum  hydrogenic  orbitals 
for  the  two  atoms,  here  designated  as  a  and  j3.  with  R  —  \Ra  —  Rp j  the  interatomic 
separation.  The  basis  of  Eq.  (70)  spans  both  spatially  symmetric  and  antisymmetric 
representations  of  the  group  S2  once  and  only  once,  whereas  its  antisymmetric  form  [Eqs 
(34)  and  (35)]  will  be  2-fold  redundant,  as  is  seen  from  Eq.  (36).  To  avoid  dealing  explicitly 
with  continuum  hydrogenic  states,  denumerable  representational  basis  sets  are  employed 
in  the  calculations  in  the  usual  way.39,40  Specifically,  approximations  to  the  orbital  spectra 
<£(“)(  1  :  Ra)  and  <p(P\ 2  :  Rp)  are  devised  which  accurately  reproduce  the  lowest-lying 
orbitals  and  which  also  span  the  Rydberg  and  lower-energy  portion  of  the  continuous 
spectrum  in  the  form  of  pseudostates.70  By  adopting  even-tempered  Gaussian  basis  sets 
for  this  purpose,82  standard  computational  methods  can  be  employed  in  evaluating  the 
matrix  elements  required  in  forming  the  spectral-product  Hamiltonian  matrix  [Eqs.  (20) 
to  (25)]  and  in  performing  its  unitary  transformation  [Eqs.  (43)  and  (60)]. 39,40 

In  Table  1  is  shown  the  spectrum  of  atomic  energies  obtained  for  s.  p.  d  and  /  orbitals 
constructed  in  even-tempered  Gaussian  basis  sets,82  with  orbital  exponents  chosen  so  that 
the  discrete  and  lower  continuum  states  are  spanned  by  the  numbers  of  orbitals  shown  for 
each  angular  momentum  value.83-85  Although  systematic  studies  of  orbital  selection  are 
not  reported  here,  it  should  be  noted  that  considerable  experience  has  been  gained  in  this 
connection  through  previous  studies  of  square-integrable  representations  of  all  the  discrete 
and  continuum  states  of  atoms  and  polyatomic  molecules  constructed  on  a  common  basis.70 
The  spectra  of  Table  1  are  judged  on  basis  of  their  spatial  and  spectral  extent  to  be  suitable 
for  describing  the  charge  distortions  accompanying  chemical  bond  or  antibond  formation  in 
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H2  in  the  interval  R  ~  1  to  5  oo,  for  describing  wave  function  symmetry  or  antisymmetry  in 
the  absence  of  explicit  electron  exchange  terms,  and  for  otherwise  approximating  spectral 
closure  in  this  interval. 

In  Figure  1  are  shown  selected  eigenvalues  of  the  (1,296  X  1,296)  metric  matrix  S  (R) 
of  Eq.  (39)  for  H2  evaluated  employing  only  the  s  and  p  basis  sets  of  Table  1. 84,85  The 
normalized  two-electron  antisymmetrizer, 

P U(i,  2)  =  (l/2)-1^2(l  —  P i2),  (71) 

is  employed  in  this  case,  with  the  redundancy  factor  Q=(2!/l!l!)=2  indicating  that  the 
totally  antisymmetric  representation  of  is  spanned  twice  by  the  antisymmetric  spectral- 
product  basis  in  the  limit  of  closure.  Only  the  fifty  largest  (sj  ~  2)  and  the  fifty  smallest 
( Si  ~  0)  eigenvalues  of  S(  R)  are  shown  in  Figure  1  as  functions  of  the  interatomic  separation 
in  H2.  The  eigenstates  of  S  (i?)  corresponding  to  Sj  ~  2  values  refer  to  linearly-independent 
spatially  antisymmetric  states  constructed  in  the  antisymmetric  basis,  whereas  those  cor¬ 
responding  to  Si  ~  0  values  refer  to  linearly-dependent  combinations  of  the  antisymmetric 
basis,  in  accordance  with  Eqs.  (43)  to  (48).  On  the  other  hand,  the  states  obtained  in  the 
spectral-product  basis  corresponding  to  the  Sj  ~  2  of  values  of  S (i?)  refer  to  antisymmetric 
states,  while  those  corresponding  to  Si  ~  0  values  refer  to  symmetric  states,  in  accordance 
with  Eqs.  (53)  to  (58).  When  the  two-electron  symmetric  projector  is  employed  in  place 
of  Eq.  (71),  results  identical  to  the  foregoing  are  obtained,  but  with  the  states  previously 
corresponding  to  Sj  ~  2  and  st  ~  0  interchanging  their  identities 

A  significant  number  of  the  eigenvalues  of  S  (R)  depicted  in  Figure  1  evidently  maintain 
their  extreme  values  (sj  ~  0  or  2)  over  the  chemical  interaction  region  (R  m  1  to  5  ao), 
whereas  only  a  very  few  of  these  extend  into  the  van  der  Waals  region  (i?  «  5  to  10  ao). 
These  behaviors  are  entirely  in  accord  with  the  spatial  characteristics  of  the  spectral  states 
of  Table  1  employed  in  constructing  S(R).  which  states  have  relatively  small  amplitudes 
at  distances  ~  5  to  10  ao  from  the  atomic  origins.  Accordingly,  the  basis  of  Table  1  can  be 
expected  to  give  converged  results  in  the  chemical  interaction  region,  and  particularly  at 
the  equilibrium  interatomic  separation  (i?  =  1.40  ao),  whereas,  alternative  or  larger  basis 
sets  are  required  to  achieve  closure  at  larger  interatomic  separations. 

The  spectra  of  Table  1  are  employed  in  calculations  of  energies  and  expectation  values 
for  the  lowest-lying  1E+  and  3E+  states  in  H2  at  the  equilibrium  interatomic  separation 
following  the  developments  of  Sections  4  and  5. 84,85  These  values  for  H2  can  be  obtained 
directly  from  diagonalization  of  the  full  spectral-product  Hamiltonian  of  Eq.  (20),  or 
from  the  unitary  transformation  of  Eq.  (33),  with  identical  values  resulting  from  the  two 
procedures  in  the  limit  of  closure.  The  total  energies  and  binding  energies  and  the  ex¬ 
pectation  values  of  the  electron  transposition  operator  Pi2  for  both  states  shown  in  Table 
2  evidently  converge  monotonically  in  both  cases  to  known  values86  with  increase  in  the 
angular  momentum  character  of  the  atomic  basis  states  employed.  Similarly,  the  norms 
of  the  familiar  singlet  and  triplet  Heitler-London  functions87  represented  in  the  spectral- 
product  basis  also  shown  in  Table  2  give  additional  indication  of  the  closure  achieved  for 
exchange  terms  in  the  wave  functions  this  case.  It  is  found  in  the  finite  basis  sets  of  Table 
1  that  significantly  more  than  one-half  of  the  transformed  Hamiltonian  matrix  of  Eq.  (33) 
is  required  in  order  to  reproduce  accurately  the  results  obtained  from  the  full  spectral- 
product  Hamiltonian  matrix.  That  is,  the  symmetric  and  antisymmetric  subspaces  in  this 
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case  are  not  completely  isolated  by  the  unitary  transformation  of  Eq.  (33)  into  approx¬ 
imately  one-half  of  the  dimension  of  the  original  spectral-product  basis,  and  the  blocks 
of  the  Hamiltonian  matrix  of  Eq.  (60)  identified  as  {Ug(i2)  •  H(i2)  •  ^Js(R)}PP/uu  must 
include  a  sufficient  number  of  small  but  non-zero  “off-diagonal”  contributions  in  order  to 
achieve  eigenvalues  in  accord  with  those  obtained  from  the  full  Hamiltonian  matrix.  These 
results  serve  to  emphasize  that  the  development  of  Section  5  provides  a  formal  proof  of 
convergence  of  the  atomic  spectral  method,  rather  than  an  optimal  computational  imple¬ 
mentation  of  the  approach,  and  that  additional  Gaussian  or  alternative  basis  functions  and 
other  computational  strategies  are  required  to  improve  the  rate  and  degree  of  convergence 
achieved. 

6.2  Recursive  Construction  of  the  Subspace  Hamiltonian 

The  metric  matrix  and  the  corresponding  explicit  unitary  transformation  of  the  Hamil¬ 
tonian  of  Section  5  are  avoided  employing  a  recursive  method  to  generate  sequences  of 
(Krylov-Lanczos)  basis  states  which  span  the  two  permutation-symmetry  subspaces  of 
the  spectral- product  basis  of  Eq.  (70).  The  development  employs  symmetric  or  anti- 
symetric  starting  functions  in  its  implementation,  essentially  constructing  the  subspace 
(i;  2  :  R )  of  Eq.  (55)  contained  in  Eq.  (70)  around  a  chosen  function  of  good  symme¬ 
try,  and  avoiding  the  unwanted  subspace  3>g^(i;  2  :  R)  entirely.  In  this  way,  tri-diagonal 
Hamiltonian  matrices  of  much  lower  dimensions  than  that  of  Eq.  (70)  can  be  constructed 
which  converge  to  unitary  equivalent  representations  of  the  totally  symmetric  or  antisym¬ 
metric  blocks  of  H(jR)  without  explicit  evaluation  and  storage  of  the  entire  metric  matrix 
or  the  Hamiltonian  matrix  required  in  Section  5. 

Following  largely  standard  procedures,80’81  an  orthonormal  symmetric  or  antisymmet¬ 
ric  Krylov-Lanczos  [vj{ i;  2  :  R ),  j  =  1,2,..  .jmax]  subspace  of  the  spectral-product  basis 
is  written  in  the  form 


Vj( i;a  :  R)  =  $(i;2  :  R)  -Vj{R),j  =  1,  2, . . .  jmax,  (72) 

where  the  column  vectors  Vj(R )  are  obtained  from  the  recursive  Lanczos  equations80’81 

Pj(R)  Vj+i(R)  =  [H(i?)  -  ctj(R)  I  ]  •  Vj(R)  -  Pj-i(R)  vj-1(R)i  j  =  1,  2, . .  .jmax-  (73) 

Here,  v0(R )  can  be  taken  as  a  null  vector  without  loss  of  generality  and  tq (R)  obtained 
using  Eq.  (72)  from  a  suitably  chosen  symmetric  or  antisymmetric  starting  function  (see 
below).  The  coefficients  oij(R)  and  (3j(R),  given  by  the  expressions 

(74) 

=  (75) 

obtained  directly  from  Eq.  (73)  and  the  imposed  orthonormality  of  the  Vj(R),  provide 
the  diagonal  and  off-diagonal  elements  of  the  symmetric  or  antisymmetric  block  of  the 
Hamiltonian  matrix  in  tri-diagonal  form.  A  sequence  of  approximations  to  either  of  the 
subspace  Hamiltonian  matrices  is  obtained  in  this  way  with  increasing  order  jmaa;,  whereas 
convergence  is  achieved  in  the  limit  jmax  oo.70 
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To  ensure  that  the  sequence  of  basis  states  of  Eqs.  (72)  to  (75)  generates  a  proper 
subspace  of  the  spectral-product  basis,  and  that  the  tri-diagonal  Hamiltonian  matrix  of 
Eqs.  (74)  and  (75)  so  constructed  is  unitarily  equivalent  to  the  symmetric  or  antisymmet¬ 
ric  block  of  the  spectral-product  Hamiltonian  matrix,  the  function  required  to  start  the 
recurrence  of  Eqs.  (73)  to  (75)  must  be  explicitly  symmetric  or  antisymmetric  and  have 
nonzero  projection  on  the  entire  spectral-product  basis.  These  requirements  ensure  that 
all  the  states  in  the  sequence  of  Eq.  (72)  will  be  explicitly  symmetric  or  antisymmetric  in 
the  closure  limit,  consequent  of  the  permutation-symmetric  nature  of  the  Hamiltonian  in 
Eq.  (73). 

Once  the  Krylov-Lanczos  sequence  has  been  constructed  employing  Eqs.  (72)  to  (75), 
corresponding  Schrodinger  eigenfunctions  are  obtained  in  the  form70 

Jmax 

Tr(i;  a  :  R)  =  NV{R)  ^  qj(Er(R))vj( i;  2  :  R ) 

3= 1 

Jmax 

=  £  qj(ET(R))Vj(R),  (76) 

3  =  1 

where  Nr(R)  =  {Yl'jZT  Qji^riR))2}-1^2  is  a  normalization  factor,  the  orthogonal  poly¬ 
nomials  qj(E )  satisfy  the  recurrence  relation 

Pj(R)  Qj+ i(E)  =  [E  -  olj  (R)  ]qj(E)  -  Pj-i(R)  qj-i(E),  j  =  1,2,..  .jmax,  (77) 

and  the  starting  conditions  for  the  recurrence  are  qo(E)  =  0  and  q\{E )  =  1.  The  associated 
energies  Er{R)  are  obtained  from  the  secular  equation  qjmax+1(Er(R ))  =  0,  with  the 
qj(E)  characteristic  polynomials  of  the  Hamiltonian  in  the  Krylov-Lanczos  basis.80,81  The 
development  of  Eqs.  (73)  to  (77)  is  identical  to  use  of  the  unitary  transformation  of  Eqs. 
(53)  to  (62)  when  jmax  approaches  the  dimension  of  the  spectral-product  basis,  but  is 
much  more  rapidly  convergent  to  specific  portions  of  the  subspace  of  physical  solutions 
than  is  the  approach  described  in  Section  5.  This  is  because  appropriately  chosen  starting 
functions  tq(i;2  :  R)  can  act  as  points  of  spectral  concentration  about  which  optimal 
Hilbert  space  descriptions  of  the  desired  physical  states  are  built  by  the  iterations  of  Eqs. 
(73)  to  (75). 70,80,81 

In  Figure  2  are  shown  XE+  and  3E+  potential  energy  curves  in  H2  obtained  from  the 
recursive  development  of  Eqs.  (72)  to  (77)  and  the  basis  of  Table  1  employing  symmetric 
and  antisymmetric  Heitler-London  test  functions  [iq(i;2  :  R)]  to  construct  the  starting 
vector  Vi(R)  in  each  case.  These  reference  states  are  excellent  zeroth-order  approximations 
to  the  physical  states  desired,  with  the  leading  tri-diagonal  Hamiltonian  matrix  element 
(a;i)  the  Heitler-London  energy  in  each  case.  As  the  iteration  of  Eqs.  (73)  to  (75)  proceeds, 
the  relatively  low-order  tri-diagonal  matrices  obtained  provide  rapidly  converging  approx¬ 
imations  to  the  desired  subspace  energy  and  wave  function  corresponding  to  the  chosen 
starting  function,  avoiding  explicit  diagonalizations  of  large-dimensioned  Hamiltonian  or 
metrix  matrices.  Evidently,  the  potential  curves  of  Figure  2  obtained  converge  rapidly  in 
the  chemical  region  (R  =  1  to  5  ao)  with  increase  in  the  orbital  angular  momentum  values 
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employed,  the  [sp]  limit  already  providing  «  90%  of  the  chemical  bonding  energy  at  the 
equilibrium  interatomic  separation,  and  the  higher  values  providing  results  in  increasingly 
good  agreement  with  the  correct  values.  The  Heitler-London  functions  for  H2  are  well 
represented  in  the  spectral-product  basis  [Eq.  (72)]  at  the  equilibrium  internuclear  sepa¬ 
ration  [see  Table  2],  but  must  be  included  explicitly  as  supplements  to  the  basis  at  larger 
separation  in  Figure  2,  rather  than  employing  the  expansion  of  Eq.  (72).  In  this  latter 
event,  one  additional  row  and  column  is  added  to  the  matrix  H(i?)  in  Eq.  (73),  although 
the  development  of  Eqs.  (72)  to  (77)  is  otherwise  unchanged.84 

Finally,  although  the  present  implementation  of  the  development  was  not  devised  specif¬ 
ically  to  study  the  van  der  Waals  region  of  bonding,  it  is  satisfying  that  the  correct  long- 
range  energy  limit  is  nevertheless  obtained  automatically  from  the  recursive  projection 
procedure.  In  Figure  3  are  shown  energy  curves  obtained  from  the  present  development 
in  the  long-range  region  (R  >  12  uq).  in  which  interval  the  singlet  and  triplet  energies 
are  indistinguishable,  in  comparison  with  values  calculated  employing  standard  singles 
and  doubles  configuration-interaction  methods  in  the  indicated  basis  sets.  These  results 
are  plotted  in  a  manner  [R°{E(R  — >  00)  —  E(R)}  vs  1  / R2]  that  allows  extraction  of  the 
two  leading  van  der  Waals  coefficients.  The  Cq  coefficients  («  6.45  au )  obtained  from 
the  (1  /R2)  — )■  0  limits  of  both  sets  of  calculations  in  the  [6s6p]  basis  set  are  evidently  in 
mutual  agreement  and  also  in  general  accord  with  the  known  correct  value  (6.499  au).89 
Similarly,  the  Cs  coefficients  (?^  123  au)  obtained  from  the  slopes  of  the  curves  in  Figure 
3  calculated  in  the  [6s6p6d]  basis  set  are  also  in  accord  with  the  known  value  (124.4  au).90 
By  contrast,  although  the  energy  curves  obtained  from  the  atomic  spectral  development 
follow  the  conventional  singles  and  doubles  configuration-interaction  calculations  for  larger 
R  (>  15  ao)  values,  they  do  not  reproduce  the  curvatures  in  the  configuration-interaction 
results  for  smaller  R  values.  This  circumstance  is  a  consequence  of  the  need  for  larger  or 
alternative  basis  sets,  and  possibly  higher  angular  momentum  terms  [see  Table  2],  required 
by  the  spectral  method  to  insure  convergence,  which  terms  are  absent  in  the  results  of  Fig¬ 
ure  3.  That  is,  the  energy  curve  in  the  12  <  R  <  15  ao  region  is  not  accurately  represented 
by  the  asymptotically  correct  Heitler-London  term  and  the  [6s6p]  and  [6s6p6d]  spectral- 
product  basis  sets  employed  in  Figure  3,  emphasizing  the  well-known  sensitivity  of  the 
small  interaction  energies  in  the  long-range  region  to  the  effects  of  electron  exchange.28’86 

7.  Discussion  and  Concluding  Remarks. 

The  atomic  spectral  method  described  here  adopts  a  new  perspective  on  calculations  of 
the  electronic  structures  of  molecules  and  other  atomic  aggregates.  Enforcement  of  wave 
function  antisymmetry  in  an  atomic  spectral-product  basis  is  deferred  until  after  construc¬ 
tion  of  the  Hamiltonian  matrix,  the  latter  consequently  taking  a  particularly  simple  form 
which  involves  only  unperturbed  atomic  energies  and  pairwise-atomic  (Coulombic)  inter¬ 
actions.  The  required  interaction  energies  are  obtained  from  atomic  transition-density 
matrices  which  can  be  calculated  once  and  for  all,  providing  thereby  a  quantitative  ap¬ 
proach  to  the  quantum  theory  of  matter  in  terms  of  spectral  attributes  of  its  atomic 
constituents.  Antisymmetry  requirements  are  subsequently  enforced  employing  a  unitary 
transformation  implemented  by  recursive  methods  devised  specifically  for  this  purpose. 
These  and  related  apspects  of  the  method  are  illustrated  with  calculations  on  molecular 
hydrogen  (H2),  which  demonstrate  convergence  to  results  obtained  from  conventional  prior 
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antisymmetrization  methods.  Additional  applications  of  the  formalism  presently  underway 
should  help  to  demonstrate  the  viability  of  the  development  in  cases  in  which  non-Pauli 
representations  are  present. 

The  theoretical  development  adopts  the  methodologies  and  attributes  of  a  variety  of 
conventional  methods  for  determining  the  adiabatic  electronic  wave  functions  and  ener¬ 
gies  of  chemical  aggregates.  Use  of  the  atomic  spectral-product  basis  dates  from  early 
perturbation-theoretrical  treatments  of  atomic  interactions25  and  related  later  efforts  to 
incorporate  the  effects  of  electron  exchange  in  this  context.26-38  In  contrast  to  these  ap¬ 
proaches,  the  present  development  is  non-perturbative  in  nature  and  is  applicable  to  both 
chemical  and  physical  (noncovalent)  interatomic  interactions  on  a  common  basis.  Although 
explicit  electron  exchange  is  avoided  at  the  outset,  the  totally  antisymmetric  subspace  of 
the  atomic  spectral-product  basis  is  isolated  and  the  unphysical  non-Pauli  representa¬ 
tions  are  avoided.26-38  Convergence  is  ensured  in  this  way  to  the  physically  significant 
eigenstates  contained  in  the  spectral-product  representation,  providing  non-perturbative 
methods  for  avoiding  the  non-Pauli  states  that  accompany  the  change  in  symmetric  group 
consequent  of  aggregate  formation.32’37  Although  convergence  is  ensured  in  the  limit  of  clo¬ 
sure,  it  remains  to  illustrate  in  greater  detail  the  computational  efficiency  of  the  method 
described  for  recursive  isolation  of  the  physical  subspace  in  more  complex  cases  in  which 
large  numbers  of  unphysical  discrete  and  continuum  states  are  present.  Additionally, 
although  no  explicit  mention  is  made  here  of  perturbative  solutions  in  the  totally  antisym¬ 
metric  subspace,  this  option  has  been  previously  considered  and  remains  to  be  explored  in 
complete  detail.45  Of  particular  interest  in  this  connection  is  the  efficient  incorporation  of 
charge-transfer  and  other  explicitly  antisymmetric  terms  in  the  form  of  chemically  signif¬ 
icant  reference  functions  to  accelerate  convergence  of  the  development  to  antisymmetric 
states,  and  to  ensure  that  the  correct  numbers  of  degenerate  multiplets  are  obtained  in 
long-range  interaction  limits. 

An  explicitly  antisymmetrized  atomic  spectral-product  basis  has  been  employed  pre¬ 
viously  in  so-called  atoms-in-molecules  approaches,46-50  dating  from  early  observations 
of  the  inadequacy  of  atomic-orbital  representations  in  describing  atomic  multiplet  order¬ 
ings  and  of  their  consequent  unsuitability  for  calculations  of  total  molecular  energies.46  The 
present  development  generalizes  the  notions  of  this  approach  in  a  rigorous  manner  employ¬ 
ing  the  atomic  charge  and  transition  densities  of  the  non-iteracting  atomic  constituents. 
As  a  consequence  of  Brillouin’s  Theorem,63  such  one-electron  properties  are  commonly 
understood  to  be  less  sensitive  to  the  effects  of  electron  correlation  than  are  total  molec¬ 
ular  energies.  Once  the  atomic  charge-  and  transititon-density  matrices  are  determined 
employing  conventional  atomic-structure  methods,  they  can  be  employed  repeatedly  in 
calculations  of  the  required  (Coulombic)  atomic  interactions.  The  incorporation  of  experi¬ 
mental  atomic  energies  in  the  Hamiltonian  matrix,  if  required,  is  straightforward,  and  the 
form  of  the  matrix  obtained  is  considerably  simpler  than  in  the  original  atoms-in-molecules 
approaches.46-50  Most  important,  the  convergence  of  the  present  method  is  assured  in 
the  spectral  closure  limit,  and  difficulties  associated  with  potential  linear-dependence  are 
avoided  entirely  in  the  orthonormal  representation  employed.  In  light  of  these  attributes, 
largely  consequent  of  separation  of  Hamiltonian  matrix  construction  from  the  enforcement 
of  wave  function  antisymmetry,  the  present  development  provides  an  atoms-in-molecules 
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strategy  alternative  to  that  devised  many  years  ago  and  employed  in  various  forms  in  the 
interim.  Although  no  explicit  mention  is  made  here  of  incorporating  diatomic  calculations 
or  other  related  information  in  aggregate  studies,51-58  this  approach  has  been  described 
previously  and  employed  in  computational  applications  to  the  structures,  spectra,  and 
dissociation  dynamics  of  atomic  clusters.45’91  The  important  theorems  upon  which  these 
studies  are  based  and  a  more  detailed  accounting  of  this  variant  of  the  atomic  spectral 
development  are  subjects  of  separate  reports. 

The  relatively  simple  forms  of  the  two-center  interaction  terms  required  in  construct¬ 
ing  the  atomic  spectral-product  Hamiltonian  matrix  allows  its  representation  in  a  wider 
range  of  orbital  functions  than  the  Gaussians  commonly  employed  in  conventional  elec¬ 
tronic  structure  calculations.  In  particular,  Slater  forms  allow  for  efficient  evaluations  of 
the  required  interaction  integrals,71,72  potentially  providing  better  representations  of  the 
atomic  states  and  associated  one-electron  transition  densities  required  in  the  development 
than  those  obtained  in  Gaussian  representations.  The  ab  initio  methods  currently  in  use 
for  constructing  the  ground  and  excited  electronic  states  of  atoms  and  molecules  can  be 
adopted  for  this  purpose  with  little  modification.7-12  Although  these  methods  as  currently 
employed  can  provide  highly  accurate  results  for  atomic  systems,  they  require  allocations  of 
significant  computational  resources  in  order  to  perform  the  repeated  electronic  structure 
calculations  at  different  molecular  geometries  required  in  constructing  potential  energy 
surfaces  and  related  properties  of  molecules.  Additionally,  they  generally  deal  with  total 
electronic  energies,  rather  than  with  atomic  binding  energies,  in  spite  of  the  well-known 
difficulties  inherent  in  such  approaches,  the  constituent  atoms  making  their  presence  appar¬ 
ent  only  through  the  commonality  of  atomic  basis  sets  employed  in  such  calculations.39-44 
By  contrast,  the  present  approach  avoids  repeated  high-level  molecular  calculations  of  to¬ 
tal  aggregate  energies  and  adopts  the  ab  initio  methods  only  in  accurate  calculations  of 
atomic  spectral  states  and  their  one-electron  transition  densities.  Use  of  the  orthonormal 
spectral-product  basis  also  circumvents  the  well-known  basis-set  superposition  difficulties 
associated  with  conventional  approaches,92  which  have  proved  particularly  troublesome  in 
studies  of  non-covalent  interactions.  Of  course,  systematic  methods  must  be  devised  for 
construction  and  tabulation  of  the  atomic  information  required  for  implementation  of  the 
present  approach,  in  which  connection  it  should  be  noted  that  atomic  spectra  are  of  con¬ 
siderable  interest  in  their  own  right,93  separate  from  their  use  in  the  context  of  calculations 
of  aggregate  electronic  structures. 

The  widely  employed  density- functional-theory  (DFT)  approaches  provide  useful  struc¬ 
tural  and  other  ground-state  information  for  large  molecules  and  atomic  aggregates,13-17 
and  can  be  extended  to  studies  of  dynamical  response  functions  for  calculations  of  excited- 
state  properties.18  Although  the  present  approach  is  based  on  largely  pedestrian  theoretical 
concepts,  and  entails  no  approximations  other  than  those  of  a  computational  nature,  it 
nevertheless  has  some  theoretical  commonality  with  the  DFT  methods  in  the  shared  at¬ 
tempt  to  circumvent  complexities  associated  with  antisymmetry  requirements,  although  by 
different  means.  Specifically,  conventional  antisymmetric  methods  are  employed  in  calcu¬ 
lations  of  the  required  atomic  spectral  information,  but  the  Coulombic  atomic  interaction 
potentials  employed  in  the  present  approach  to  aggregate  formation  are  constucted  in  the 
absence  of  inter-atomic  exchange  effects,  and  are,  therefore,  largely  unphysical.  Relatedly, 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


attempts  to  incorporate  explicitly  the  subtle  effects  of  electron  correlation  in  the  DFT 
approaches  to  van  der  Waals  interactions  apparently94  lead  naturally  to  the  one-electron 
transition-density  matrices  of  particular  computational  focus  in  the  present  development, 
suggesting  that  study  of  possible  adoption  of  DFT  methods  in  the  present  context  may  be 
beneficial. 

Finally,  although  it  is  premature  to  discuss  the  possible  role  of  the  present  develop¬ 
ment  in  connection  with  semi-empirical  and  QM/MM  approaches  to  larger  systems,19-24 
there  may  ultimatly  be  some  relevance  in  this  connection  in  view  of  the  atomic-interaction 
perspective  adopted,  and  in  light  of  the  growing  importance  of  such  methods  to  studies  of 
proteins  and  other  biological  macromolecules.24  Particularly  needed  are  higher- level  QM 
methods  to  treat  proton  and  electron  transfer,  and  the  isomerizations  and  other  motions 
of  optical  chromophores,  in  the  presence  of  solvating  protein  or  other  environments.  Ad¬ 
ditionally,  it  would  be  highly  desirable  to  develop  methods  that  can  avoid  the  QM-MM 
boundary  effects  experienced  in  some  applications  of  current  methods,  among  the  other 
challenges  such  studies  present. 

The  foregoing  issues  and  additional  matters  associated  with  more  general  implementa¬ 
tions  of  the  atomic  spectral  approach  clearly  require  studies  well  beyond  the  scope  of  the 
present  report.  Investigations  of  these  matters  and  the  results  of  computational  applica¬ 
tions  to  other  systems  currently  in  progress  will  be  described  in  due  course. 
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Table  1.  Spectral  Energies  for  Atomic  Hydrogen.® 


s-basisb 

p-basisb 

d-basisb 

f-basisb 

-0.499991 

-0.124998 

-0.055256 

-0.026707 

-0.124994 

-0.053439 

-0.011199 

0.028807 

-0.048465 

0.025060 

0.108901 

0.187597 

0.088668 

0.274825 

0.450236 

0.617389 

0.563368 

0.986192 

1.361055 

1.748440 

1.951139 

5.660465 

15.152536 

38.982050 

98.333482 

249.224098 

657.488778 

2.880038 

7.854649 

21.314359 

3.708833 

9.719458 

25.786063 

4.791877 

°  Orbital  energies  (ait)  obtained  from  diagonalization  of  the  atomic 
hydrogen  Hamiltonian  employing  the  indicated  basis  sets.82-85 
b  Basis  sets  employed  are  the  most  diffuse  [12s8p8d6f]  primitives 
chosen  from  a  set  of  14  even-tempered  Gaussian  primitives  for 
each  angular  momentum  symmetry,  which  set  consists  of  12  reg¬ 
ularized  even-tempered  primitives82  supplemented  with  two  addi¬ 
tional  primitives  in  each  symmetry  having  exponents  of  0.02786 
and  0.01156. 
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Table  2.  Electron  Pair-Bond  Calculations.® 


Basis6 

Energy(a.u.)c 

Binding(eV)c 

<A2)c 

(  $(#L)  2)c 

1E+  state 

M 

-1.0096 

+0.2618 

+0.5255 

+0.8052 

M 

-1.0691 

+1.8809 

+0.8409 

+0.9469 

[spd] 

-1.1140 

+3.1027 

+0.9525 

+0.9872 

[spdf] 

-1.1384 

+3.7667 

+0.9847 

+0.9967 

Exact86 

-1.1745 

+4.7478 

+1.0000 

+1.0000 

3E+  state 

M 

-0.5586 

-12.0109 

-0.6226 

+0.2950 

M 

-0.6641 

-9.1400 

-0.6905 

+0.8078 

[spd] 

-0.7249 

-7.4856 

-0.9317 

+0.9538 

[spdf] 

-0.7524 

-6.7372 

-0.9801 

+0.9870 

Exact86 

-0.7842 

-5.8737 

-1.0000 

+1.0000 

a  Values  at  R  =  1.40  «o  obtained  from  diagonalization  of  the 
Hamiltonian  matrix  of  Eqs.  (20)  to  (25)  for  H2,  or,  equivalently, 
from  the  unitary  transformation  of  Eq.  (33)  in  the  text.84,85 
b  Denotes  the  portion  of  the  [12s8p8d6/]  basis  set  indicated  in 
Table  1  employed  in  the  calculation. 

c  Total  and  binding  energies  as  indicated;  (Pi2)  refers  to  the  expec¬ 
tation  value  of  the  electron  transposition  operator  i\2;  |<1>(-HX)|2 
is  the  norm  of  the  Heitler-London  function  as  represented  in  the 
spectral-product  basis. 
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Captions  for  Figures 


Figure  1.  Eigenvalues  s*  of  the  metric  matrix  of  Eqs.  (39)  and  (71)  for  H2,  constructed 
in  the  [12s8p]  basis  set  of  Table  1  as  functions  of  interatomic  separation  i?(oo):84,85  val¬ 
ues  Si  ~  2  refer  to  approximately  antisymmetric  eigenfunctions  of  electron  coordinates 
constructed  in  the  spectral-product  basis,  whereas  values  s*  ~  0  refer  to  approximately 
symmetric  functions  of  electron  coordinates,  as  is  discussed  more  fully  in  the  text;  apparent 
discontinuities  in  the  slopes  of  the  plots  are  due  to  the  coarse  grid  in  R  values  employed. 

Figure  2.  Total  energy  ( au )  curves  for  the  1E+  and  3E+  states  of  H2  as  functions  of  atomic 
separation  f?(ao):  heavy  solid  lines  refer  to  previously  determined  accurate  values;86  light 
solid  lines  refer  to  Heitler-London  values;87  dashed  lines  give  the  present  results  obtained 
from  the  recursion  procedure  indicated  in  the  text  employing  the  [s],  [sp],  and  [spd\  basis 
states  of  Table  1  and  Heitler-London  test  functions  in  each  case.83,84 

Figure  3.  Potential  energy  ( au )  curves  for  the  1E+  state  of  molecular  hydrogen  in  the  long- 
range  interaction  region  obtained  as  indicated:  open  squares  -  accurate  values  obtained 
from  single  and  doubles  configuration-interaction  calculations  in  [6s6p]  and  [6s6p6d]  basis 
sets  consisting  of  six  regularized  even-tempered  primative  Gaussian  orbitals  of  each  angular 
momentum  symmetry;82  solid  circles  -  spectral  theory  results  obtained  from  the  recursive 
development  of  Eqs.  (73)  to  (77),  calculated  employing  the  indicated  basis  sets;  lines  - 
fits  to  the  straight-line  portions  of  the  calculated  values;  the  (1/R2)  — *  0  intercept  gives 
Ce  —  6.45  au,  whereas  the  slope  of  the  [6s6p6d]  curve  in  this  limit  gives  Cg  =  123  au,  in 
good  accord  with  the  known  values  of  6.499  and  124.4  au,  respectively.89,90 
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Eigenvalue 


[EH-E(R)]  (au) 
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